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EXECUTIVE  SUMMARY 


The  increasing  demand  for  multimedia  interactive  services  and  great  mobility  led 
the  various  communications  industries  to  develop  the  third  generation  of  service  (3G). 
While  the  existing  second-generation  system  was  originally  designed  for  voice 
transmissions  only,  the  third  generation  supports  numerous  high  data  rate  applications 
such  as  full  wireless  Internet  access.  The  second-generation  system  can  reach  up  to  10 
kbps  while  the  third  generation  can  reach  up  to  2  mbps.  In  addition,  the  third  generation 
wireless  system  guarantees  144  kbps  for  faster  moving  users.  The  existing  systems  use 
different  multiple  techniques  like  Time  Division  Multiple  Access  (TDMA)  and 
Frequency  Division  Multiple  Access  (FDMA),  which  are  bandwidth  limited.  The  next 
generation  of  wireless  systems  (3G)  uses  Wideband  Code  Division  Multiple  Access  (W- 
CDMA).  Thus,  the  system  capacity  is  increased  many  times  compared  to  the  existing 
systems.  This  type  of  multiple  access  is  interference  limited.  For  this  reason,  W-CDMA 
systems  have  to  use  different  types  of  techniques  to  overcome  intracell  and  intercell 
interference  as  well  as  the  additive  white  Gaussian  noise  (AWGN). 

Accordingly,  a  forward  channel  for  a  DS-CDMA  cellular  system  is  created.  The 
information  signal  is  developed  which  propagates  through  the  communication  channel. 
The  extended  Hata  model  is  being  used  in  order  to  compute  the  large-scale  path  loss. 
Furthermore,  since  the  transmitting  signal  propagates  through  different  types  of  terrains, 
the  Nakagami  random  variable  is  introduced  in  order  to  incorporate  small-scale  fading 
conditions.  Fluctuations  to  the  received  power  at  points,  which  have  the  same  distance 
from  the  transmitter,  are  considered  also  by  the  use  of  Lognonnal  Shadowing.  Taking 
into  account  all  these  propagation  effects  to  the  information  signal,  the  result  is  the 
received  signal  to  the  intended  mobile  user’s  receiver. 

First,  the  intended  user  is  placed  at  the  corner  of  the  hexagon  of  the  center  cell  of 
a  seven-cell  structure,  which  is  the  worst-case  scenario.  Given  an  upper  bound  on  the 
probability  of  bit  error  using  Forward  Error  Correction  (FEC),  the  performance  of  the 
cellular  system  is  analyzed.  A  statistical  model  is  created  in  order  to  approximate  the 


sum  of  d  multiplicative  Nakagami-square-Lognormal  random  variables  as  a 
multiplicative  Nakagami-square-Lognormal  random  variable.  Our  model  is  being  tested 
in  different  fading  and  shadowing  conditions.  Modeled  results  are  compared  with  those 
predicted  by  Monte  Carlo  simulations. 

In  order  to  implement  more  realistic  and  accurate  results,  the  user’s  distribution  is 
incorporated  in  the  cell  rather  than  placing  the  user  at  the  corner  of  the  hexagon  which  is 
the  worst-case  scenario.  Furthermore,  fast  power  control  is  incorporated  in  order  for  the 
received  power  by  the  mobile  user  to  be  at  a  target  level  and  thus  eliminating  the  effect  of 
Lognormal  Shadowing. 
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I.  INTRODUCTION 


A.  BACKGROUND 


The  increasing  demand  for  multimedia  and  interactive  services  led  the  various 
communications  industries  to  develop  the  third  generation  of  service  (3G).  The  existing 
second-generation  system  is  widely  being  used  mainly  for  voice  transmissions  supporting 
data  rates  up  to  10  kbps.  However,  many  applications  like  full  wireless  Internet  access 
demand  much  higher  data  rates.  The  third  generation  of  service  can  achieve  data  rates  up 
to  2  mbps  and  for  faster  users,  up  to  144  kbps.  This  new  system  employs  Code  Division 
Multiple  Access  (CDMA)  and  thus  is  able  to  increase  the  system  capacity  many  times 
compared  to  the  present  systems.  Japan  is  the  first  country  in  the  world,  which  already 
uses  this  new  type  of  wireless  communications.  Europe  and  the  United  States  is  expected 
to  implement  this  new  technique  in  2002,  and  2004  respectively. 


B.  OBJECTIVE 


Code  division  multiple  access  systems  are  interference  limited.  On  the  other 
hand,  Frequency  Division  Multiple  Access  (FDMA)  and  Time  Division  Multiple  Access 
(TDMA)  systems  are  bandwidth  limited  systems.  While  the  transmitted  signal  from  the 
base  station  propagates  through  our  communication  channel,  it  is  being  largely  affected 
by  interference  as  well  as  by  the  Additive  White  Gaussian  Noise  (AWGN).  In  addition, 
and  since  the  signal  propagates  through  different  types  of  terrain,  it  is  undergoing  a  slow 
flat  fading  type.  A  more  general  type  of  fading  is  Nakagami  fading.  By  varying  the 
Nakagami-/;/  parameter,  several  types  of  fading  such  as  one-sided  Gaussian,  Rayleigh 
and  Ricean  can  be  implemented.  However,  different  points  located  at  the  same  distance 
from  the  transmitting  base  station  have  vastly  different  path  loss.  This  phenomenon  is 
called  Fognormal  Shadowing.  Taking  into  account  all  the  aforementioned  facts,  the 
performance  of  a  direct  sequence  CDMA  six-sectored  cellular  system  in  a  Nakagami 
fading  and  Fognormal  Shadowing  environment  is  analyzed  by  computing  an  upper  bound 
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on  the  probability  of  bit  error  with  the  use  of  a  statistical  model  and  Monte  Carlo 
simulations.  Additionally,  by  complementing  different  techniques  such  as  user 
distribution  and  fast  power  control,  the  performance  of  the  system  is  enhanced. 


C.  RELATED  WORK 

A  lot  of  research  on  the  DS-CDMA  channel  is  on  the  reverse  link,  which  is  in 
general,  different  from  the  forward  link.  A  very  analytical  work  of  a  DS-CDMA  forward 
channel  has  been  developed  in  [1].  The  analysis  in  [1]  is  concerned  with  Rayleigh  fading 
instead  of  Nakagami  fading.  However,  many  of  the  concepts  are  the  same,  and  are  very 
useful  for  the  analysis  in  this  thesis.  Furthermore,  [2]  optimizes  power  using  pilot  tone 
power  control  in  a  Rayleigh-Lognormal  forward  channel. 

In  summary,  it  can  be  concluded  that  a  more  general  analysis  of  Nakagami  fading 
including  and  extending  previous  research  needs  to  be  accomplished. 


D.  THESIS  OUTLINE 


In  Chapter  II,  a  forward  channel  for  the  DS-CDMA  cellular  system  is  created.  By 
taking  into  account  all  the  appropriate  losses,  interferences  and  phenomena  such  as 
fading  and  Lognormal  Shadowing,  the  result  was  a  closed  formula  for  the  received  signal 
to  the  intended  user. 

In  Chapter  III,  the  intended  user  was  located  in  the  least  optimum  position  within 
the  center  cell,  and  given  the  expression  for  the  upper  bound  on  the  probability  of  bit 
error  with  the  use  of  Forward  Error  Correction  (FEC)  in  the  form  of  convolutional 
encoding  with  soft-decision  decoding,  the  perfonnance  of  the  cellular  system  is  analyzed. 
This  analysis  creates  a  statistical  model  in  order  to  approximate  the  sum  of  d 
multiplicative  Nakagami-square -Lognormal  random  variables  as  a  multiplicative 
Nakagami-square-Lognormal  random  variable.  The  model  is  being  tested  in  various 
fading  and  shadowing  environments  and  the  Nakagami-zn  parameter  is  varied.  The 
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modeled  results  are  compared  with  the  results  predicted  by  95,000  Monte  Carlo  trials. 
The  values  of  the  Nakagami-/;/  parameter  used  are:  m  =  0.5,  m  =  0.75,  in  =  1,  m  =  1.5  and 
m  =  2. 

In  Chapter  IV,  the  perfonnance  analysis  completed  in  Chapter  III  is  modified  by 
incorporating  user  distribution  within  the  cell  for  more  realistic  results.  Furthermore,  fast 
power  control  is  implemented  in  order  to  enhance  system  performance.  The  analysis  is 
done  in  different  fading  environments  by  varying  the  Nakagami-m  parameter. 

Finally,  in  Chapter  V,  the  conclusions  are  summarized  and  areas  for  further 
research  are  discussed. 
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II.  FORWARD  CHANNEL  MODEL 


The  forward  channel  in  a  Direct  Sequence  Code  Division  Multiple  Access  cellular 
system  is  the  communication  channel,  which  carries  traffic  from  the  base  station  to 
mobile  users.  The  reverse  channel  carries  traffic  from  mobile  users  to  the  base  station. 

In  this  chapter,  a  channel  model  for  a  DS-CDMA  cellular  system  including  fading 
and  shadowing  effects  will  be  built.  It  is  known  that  a  mobile  user  is  often  surrounded  by 
various  obstacles  such  as  buildings,  trees  and  mountains.  These  cause  the  presence  of 
phenomena  like  reflection,  diffraction  and  scattering  of  the  forward  transmitted  signal 
from  the  base  station. 

The  average  signal  power  that  a  mobile  user  receives  depends  on  the  distance 
between  the  transmitting  base  station  and  the  mobile  user,  as  well  as  on  the  type  of 
environment  between  them  [3].  The  transmitted  base  station  forward  signal  arrives  at  the 
mobile  user’s  receiver  via  many  paths  with  different  lengths  and  with  various  signal 
powers.  The  sum  of  the  multiple  signals,  arriving  at  the  mobile  receiver  leads  to  a  signal 
that  fluctuates  very  fast  as  the  mobile  and/or  environment  move(s).  Such  signal 
fluctuations  are  known  as  fading  [4],  The  fluctuations  that  occur  over  sub-wave  lengths 
scale  are  known  as  fast  fading  or  small-scale  fading  and  those  occurring  over  several 
wavelengths  are  known  as  slow  fading  or  large-scale  fading  [4],  Some  large  scale  and 
small-scale  propagation  models,  which  are  used  in  order  to  detennine  the  average  signal 
power  received,  will  be  discussed  later  in  this  chapter. 


A.  DS-CDMA  FORWARD  SIGNAL 


A  basic  seven-cell  cluster  as  shown  in  Figure  2.1  will  be  used  for  our  analysis. 
Each  cell  is  represented  by  a  hexagon  whose  center  is  the  position  where  the  base  station 
of  each  cell  is  located.  The  center  cell  is  the  main  interest.  Let  K  represent  the  total 
number  of  mobile  users  or  active  channels  in  the  center  cell,  while  Kj  is  the  total  number 
of  mobile  users  or  active  channels  in  the  adjacent  cells  i=  1,2. .  .6. 
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Figure  2. 1 .  Seven  -Cell  Cluster. 

S0(t)  represents  the  transmitted  signal  by  the  base  station  in  each  cell  which 

includes  the  traffic  for  all  active  channels  in  the  cell.  Due  to  the  DS  technique,  the  power 
spectrum  of  these  signals  has  been  spread  by  a  factor  of  N.  In  this  analysis,  the  type  of 
modulation  used  is  BPSK.  Suppose  that  there  is  a  user  located  at  the  center  cell  of  Figure 
2.1.  This  mobile  user  will  receive  traffic  that  is  intended  for  other  users  in  his  cell,  or  in 
our  case,  the  center  cell  of  Figure  2.1.  This  type  of  traffic  is  called  intracell  interference. 
In  addition  to  the  previous  traffic,  our  user  will  receive  signals  Sf(t)  transmitted  by  the 

base  stations  of  the  six  adjacent  cells.  This  type  of  interference  is  called  intercell 
interference,  or  co-channel  interference. 

1.  Walsh  Functions 

Walsh  functions  are  used  in  DS-CDMA  cellular  systems  in  order  to  spread  the 
traffic  over  a  finite  bandwidth  and  eliminate  intracell  interference.  In  the  forward  link, 
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the  traffic  is  multiplied  with  distinct  repeating  Walsh  sequences  that  are  assigned  to  each 
channel  for  the  duration  of  the  call.  The  Walsh  sequences  (functions)  are  mutually 
orthogonal.  In  that  way,  in  the  despreading  process,  the  mobile’s  user  receiver  is 
applying  his  own  unique  Walsh  function  and  despreads  only  the  traffic  intended  for  him. 
All  the  other  traffic  from  the  same  cell  at  the  specific  receiver  goes  to  zero.  However, 
some  intercell  interference  still  occurs  due  to  that  fact  that  multipath  and  signals  from 
other  cells  that  are  not  time-aligned  with  the  desired  signal  are  present.  To  reduce  this 
type  of  interference,  a  PN  sequence  modulation  process  is  used.  Thus,  it  can  be  stated 
that  the  multiple  access  scheme  for  the  forward  link  in  a  DS-CDMA  system  incorporated 
the  use  of  orthogonal  Walsh  sequences  (functions)  and  PN  sequences  of  different 
purposes. 

The  Walsh  functions  of  order  N  can  be  defined  as  a  set  of  N  time  functions  as 
follows:  te  (0,T),  j  =  0,1...,  A-l]  where  Wj(t)& {+1,-1}  ,  except  at  the  jumping 

points,  where  it  equals  to  zero.  The  quantity  W-(0)  =  1  for  all  j  while  Wj(t)  has  exactly  j 
sign  alternations  in  the  interval  (0,7)  [5] 


J(j  w;(twk(t)dt 


jo,  if  j*k 
|  T,  if  j  =  k 


In  Figure  2.2,  a  set  of  Walsh  functions  of  order  eight  are  shown. 
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Figure  2.2.  Walsh  Functions  of  Order  8. 
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It  is  now  possible  to  convert  the  ±1  amplitudes  of  the  Walsh  functions  to  a  binary 
logic  {0,1}  representation  with  the  following  conversions  +1  — »  "0"and-l  — »  "1"  [5]. 

Thus,  Table  2.1  is  created. 


Index 

Index 

Walsh  sequences 

integer  j 

sequence 

Of  order  8 

0 

000 

W0= 00000000 

1 

001 

Wj= 00001111 

2 

010 

W2= 00111100 

3 

011 

W3  =  001 10011 

4 

100 

W4  =  01100110 

5 

101 

W5  =  01101001 

6 

110 

W6  =  01011010 

7 

111 

W7  =  01010101 

Table  2.1.  The  Walsh  Sequences  of  Order  8  from  [5]. 

There  are  many  ways  that  Walsh  functions  can  be  generated  such  as  using 
Rademacher  functions  and  Hadamard  matrices. 

2.  Direct  Sequence  Spread  Spectrum  CDMA 

A  DS-Spread  Spectrum  technique  is  used  in  order  to  spread  the  power  spectrum 
of  a  base  band  signal  over  a  specified  transmission  bandwidth  [1].  DSSS  systems  are 
useful  in  military  applications  where  a  low  probability  of  intercept  (LPI),  low  probability 
of  detection  (LPD)  and  anti-jam  protection  are  required  [5].  In  third  generation  cellular 
systems,  DSSS  can  be  used  for  multiple  access  applications.  Given  a  finite  bandwidth, 
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each  base  station  located  in  each  cell  has  to  use  its  own  PN  spreading  signal  and  Walsh 
functions  for  orthogonal  covering.  Now  the  mobile  receiver  despreads  only  the  traffic 
intended  for  him  using  the  cell’s  PN  signal  and  his  unique  Walsh  function. 

In  order  to  generate  a  DSSS  signal,  it  is  necessary  to  multiply  the  information 
signal  by  a  spreading  signal.  Let,  as  in  [1],  bk(t)  represents  the  information  signal 

intended  for  mobile  user  k  where,  bk (t)  e  {±  1}  nT  <  t  <  (n  +  1)T  for  n  =  0,1,2...  with  a  bit 

interval  T.  Each  data  bit  is  spread  by  a  factor  of  N  using  a  spreading  signal  with  chip 
T 

duration  T  =  —  for  all  users. 

c  N 

PN  sequences  are  generated  from  linear  feedback  shift  registers.  PN  sequences 
and  Walsh  functions  are  used  to  spread  the  data  signal.  Walsh  functions  eliminate  the 
intracell  interference  while  PN  sequences,  to  a  lesser  degree,  combat  the  intercell 
interference.  As  stated  in  [1],  the  application  of  Walsh  function  W-(t) e  WN  alone  does 

not  guarantee  a  spreading  factor  of  N.  The  PN  signal  in  the  center  cell  is  represented  as 
c(t)  while  PN  signals  of  the  adjacent  cells  are  represented  as  ct (t)  for  i  =  1,2, ...6. 

Assume  that  a  user  is  in  the  center  cell.  In  order  for  this  user  to  correctly  receive  the  data 
intended  for  him,  his  receiver  chipping  signal  has  to  be  synchronized  with  the  chipping 
signal  c(t)  of  his  cell  base  station.  Since  this  mobile  user,  who  is  in  the  center  cell,  is  not 
synchronized  with  base  stations  in  adjacent  cells,  the  inter-cell  traffic  remains  spread  over 
the  bandwidth  and  minimizes  its  effect  in  SNR. 

3.  Transmitted  Signal,  S(t) 

Following  the  assumption  in  [1],  we  get: 
k  =  the  mobile  user  or  channel  k  in  the  center  cell 

bk  (t)  =  the  information  signal  consisting  of  binary  data  stream  for  the  k- th  user 
channel  in  the  center  cell 

Ptk  =  the  average  transmit  power  in  the  k- th  channel 


9 


Wk(t)  =  the  Walsh  function  of  the  k- th  user  channel  in  the  center  cell 
c(t)  =  PN  spreading  signal  for  the  center  cell 
fc  =  the  carrier  frequency  of  the  signal 

Then,  the  transmitted  signal  S0(t)  from  the  base  station  of  the  center  cell  is  given  by  [1] 
as  follows: 

*(0  =  Z  fiPkM‘Wt(t)c(t)cos(2xlt)  (2.1) 

k=0 

where  K  is  the  number  of  active  channels  in  the  center  cell. 

B.  PROPAGATION  IN  THE  MOBILE  RADIO  CHANNEL 


As  the  transmitted  signal  travels  from  the  base  station  to  the  mobile  users  it  looses 
power  as  a  function  of  the  distance  and  the  type  of  environment.  This  loss  in  signal 
power  is  called  path  loss,  which  is  represented  as  L  in  dB. 

Propagation  measurements  in  a  mobile  radio  channel  have  shown  that  the  average 
received  signal  power  Pr  at  a  distance  d  from  the  transmitter  is  given  as  follows: 


P..  =Pn 


rdyn 

\doJ 


(2.2) 


where  P0  is  a  reference  power  received  at  distance  d0  from  the  transmitter,  and  n  is  the 
path  loss  exponent  [3].  Table  2.2  shows  some  values  of  n  in  different  environments. 
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Environment 

Path  Loss  Exponent,  n 

Free  space 

2 

Urban  area  cellular  radio 

2.7  to  3.5 

Shadowed  urban  cellular  radio 

3  to  5 

In  building  line-of-sight 

1.6  to  1.8 

Obstructed  in  building 

4  to  6 

Obstructed  in  factories 

2  to  3 

Table  2.2.  Path  Loss  Exponents  in  Different  Environments  from  [3]. 


1.  Free  Space  Propagation  Model 


The  same  symbols  as  in  [1]  will  be  used: 
Pfr  =  the  free  space  power  received 

Pt  =  the  transmitted  power 
Gt  =  the  transmitter  antenna  gain 
Gr  =  the  receiver  antenna  gain 


Ac  =  the  wavelength  of  the  carrier,  x.  = 


fc 


d  =  the  separation  distance  between  transmitter  and  receiver 
Ls  =  the  system  (hardware)  loss  factor 

fc  =  the  carrier  frequency 
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Assume  for  simplicity  that  Gt=Gr=Ls=  1.  As  seen  from  Table  2.2  for  free 
space,  n  =  2.  Thus,  using  Equation  (2.2)  with  n  =  2,  it  is  then  possible  to  calculate  Pf. 


PJd)  =  PfM 


f  d,  3 


d0<d 


\u  j 


(2.3) 


where  Pfr{d0 )  is  the  free  space  reference  power  which  can  be  measured  or  calculated 
using  the  Friis  free  space  equation  given  by 


P,G,G,\ 

IAjtUI'L 


The  free  space  propagation  model  is  used  when  there  is  an  unobstructed  line-of- 
sight  path  between  the  transmitting  and  receiving  antenna  [1],  However,  usually  in  a 
mobile  environment,  many  times  there  is  no  line-of-sight  path  so  the  mobile 
communication  channel  cannot  be  characterized  by  the  use  only  of  the  free  space 
propagation  model.  The  Hata  model  is  one  of  the  other  large-scale  propagation  model, 
which  is  used  to  calculate  the  median  path  loss. 


2.  The  Hata  Model 


The  Hata  Model  is  the  analytical  expression  of  the  Okumara  empirical  model, 
which  is  one  of  the  most  accurate  empirical  models  in  cellular  communications  [3].  The 
Hata  model  computes  the  median  path  loss  in  an  urban  environment  and  supplies 
correction  equations  in  order  to  be  applicable  to  other  situations.  The  extended  Hata 
model  can  be  used  for  carrier  frequencies  from  1500  MHz  to  2000  MHz.  From  [1]  and 
[3]  it  can  be  seen  that 

Lmb  =46.3  +  33.9  log  /’  -  1 3.82  log hbase  -  a(hmobile )  +  ( 44.9  -  6.55  log hbase )  log d  (2.5) 

+C„ 
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^  ^ mobile  )  =  (1.1  \ogfc-0.l)hmobile -(1.56  log/c  -0.8)(dB) 


(2.6) 


[ 0  dB,  for  medium  sized  city  and  suburban  areas 
[3  dB,  for  metropolitan  centers 


(2.7) 


where  hhase  and  hmobile  are  the  heights  of  the  base  station  and  mobile  antennas,  fc  the 
carrier  frequency,  d  is  the  distance  between  the  base  station  and  the  mobile,  CM  is  an 
operating  area  correction  factor.  Equations  (2.5)  and  (2.6)  are  valid  when  the  following 
units  for  each  parameter  are  used:  fc  in  MHz,  hbase  and  hmobile  in  meters  and  d  is 

measured  in  kilometers  (km).  The  following  restrictions  apply  in  the  parameters,  which 
are  used  in  the  extended  Hata  Model: 


fc:  1500  MHz  to  2000  MHz 


Kase  '  30m  t0  200m 

Kobile  ■  lm  to  10m 

d:  1km  to  20km 

The  extended  Hata  Model  is  used  to  predict  the  median  path  loss  because  it  better 
matches  with  the  size  and  frequency  ranges  of  future  cellular  systems. 


3.  Lognormal  Shadowing 


The  average  path  loss  is  a  function  of  the  path  loss  exponent  n  and  is  given  in  dB 


as  follows:  LlldB(d)  =  Ltw(d0)  +  1 0«  log 


rd' 

V^o  J 


[1].  Taking  into  account  that  the  surrounding 


environment  may  differ  a  lot  at  two  separate  points,  which  have  the  same  distance  d  from 
the  transmitter,  it  can  be  stated  that  the  previous  equation  does  not  predict  with  accuracy 
the  path  loss  at  these  two  different  points  [3],  These  differences  can  be  measured  in 
decibels  and  they  can  be  represented  with  a  Gaussian  random  variable.  Thus,  the  path 
loss  with  shadowing  can  be  defined  as: 
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(2.8) 


4</S  (d)  4ffi(4  +  4ffl 


where  XdB  is  a  Gaussian  random  variable,  XdB  ~  N (0,<7//;)  with  zero  mean,  and  standard 
deviation  a dB  ,  and  LdB(d)  the  path  loss  for  distance  d. 

Since  XdB  is  a  zero  mean  Gaussian  random  variable,  LxdB  is  also  a  Gaussian 

random  variable  LxdB  ~  N  ( LdB ,  ad/j )  where  LdB  is  the  mean  and  the  median  value  since 

for  a  Gaussian  random  variable  the  mean  value  equals  the  median  value.  Therefore,  in 
our  model,  the  median  path  loss  predicted  by  the  extended  Hata  Model  will  be  used  as 
follows: 

LxdB=LHdB+XdB  (2-9) 

The  following  equations  are  from  [1]: 

4 dB  =  io  log  4 
L HdB  ~  10  lOg  4/ 


XdB  =  10  log  X 
4  =  lhx 

where  X  ~  X(Q,AadB)  is  a  Lognormal  random  variable  which  is  the  conversion  of  the 
Gaussian  random  variable  XdB  ~  N  (0,  adB ) .  The  parameters  of  X  are  ,u.  =  AjUdB  =  0  and 
(7x  =  A(TdH  where  x  =  tn '  0^^  [2],  Lx  is  also  a  Lognormal  random  variable 

4  ~  X  (^LHdB  ■>  dB  )  • 


It  is  known  that  the  received  power  is  given  by 

_  PtGtGr  _  Pt 
”  L  4  L 


(2.10) 


since  G \  =  Gr  =Ls=  \  have  been  nonnalized  to  one. 

However,  in  our  model,  it  will  sometimes  be  necessary  to  vary  the  value  of  Pt . 
Thus,  by  using  Pt  as  our  main  value  of  transmitting  power  at  most  of  our  channels,  we 
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can  define/)  k  as  the  signal  power  assigned  to  each  channel  k  using  a  power  factor  fk  as 
follows: 


PtJC  =fkp. 


(2.11) 


Taking  into  account  our  Hata-Lognonnal  model  for  path  loss,  L=LHX,  the 
received  power  in  the  k- th  channel  can  be  converted  into  a  random  variable  Pk  as 


fkP, 


Lh  (d)X 


(2.12) 


where, 


f k  =  the  power  factor  used  to  adjust  the  power  in  the  k- th  channel 


Pt  =  the  baseline  signal  power 


Lh  =  the  median  path  loss  using  the  Hata  Model 
X  =  the  Lognormal  random  variable  A  ( 0.  act dlj ) 

From  [1]  it  is  shown  that  the  inverse  of  a  Lognonnal  random  variable  is  also  a 
Lognormal  random  variable.  Thus,  the  power  received  in  the  k- th  channel  in  Equation 
(2.12)  is  a  Lognormal  random  variable,  Pk  ~  A  [jupi ,  A(7dfj )  where  flpt  =  tn(flPliL„). 

The  model  that  was  developed  in  this  Section  takes  into  account  large-scale 
propagation  losses  and  Lognormal  shadowing. 


4.  Small-Scale  Fading  due  to  Multipath 

The  signal  fluctuations  that  occur  over  sub-wavelengths  have  already  been 

defined  and  are  known  as  fast  fading  or  small-scale  fading  [4],  due  to  the  multipath 

reception  of  a  wireless  signal.  Many  signal  copies  arrive  at  the  receiver  at  different  time 

intervals  because  they  follow  different  propagation  paths  or  because  in  a  wireless 

environment  either  transmitter  or  receiver  moves  with  time.  This  is  why  there  is  a 

Doppler-shift  in  the  frequency  of  the  transmitted  signal.  The  fact  that  the  receiver  will  or 
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will  not  be  able  to  receive  correctly  each  information  data  bit  or  symbol  depends  on  the 
small-scale  fading  characteristics  of  the  channel  and  the  signal  characteristics  itself.  In 
[3],  a  classification  of  types  of  small-scale  fading  occurs  due  to  two  independent  facts 
such  as  multipath  delay  spread  and  Doppler  spread.  Figure  2.3  shows  a  tree  of  the  four 
different  types  of  fading  [3], 


Figure  2.3.  Types  of  Small-Scale  Fading  from  [3]. 

It  can  also  be  mentioned  that  in  reality,  fast  fading  environments  have  only  been 
observed  in  communication  systems  with  very  low  data  rates  [3].  Since  high  data  rates 
are  going  to  be  used,  our  model  is  characterized  as  slow  fading  [1].  A  general  method  for 
modeling  the  amplitude  variations  in  a  flat  fading  channel  is  by  assuming  the  amplitudes 
are  distributed  as  a  Nakagami  random  variable.  In  the  next  Section,  the  effects  of  large- 
scale  path  loss  and  small  scale  fading  will  be  taken  into  account  in  our  channel  model  and 
will  be  defined  as  the  Nakagami-Lognormal  channel  model. 
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c. 


NAKAGAMI-LOGNORMAL  CHANNEL  MODEL 


The  Nakagami-Lognormal  channel  model  can  be  said  to  be  a  slow-flat-Nakagami 
fading  channel  with  Lognonnal  shadowing  and  with  path  loss  which  are  calculated  by  the 
extended  Hata  Model.  If  the  same  procedure  used  in  [1]  is  followed,  and  the  Rayleigh 
random  variable  is  changed  to  a  Nakagami  random  variable,  we  find  that  in  a  DS-CDMA 
cellular  system,  the  received  signal  by  a  mobile  user  in  the  center  cell  is  given  by 


KO  =  Z  R^Abk  (t)wk  (*M0  cos(2  nfct) 

6  £  (2-13) 
+1  Z  M  IPjbij  ( t  +  r,  )Wy(t  +  rt ) c,  (; t  +  t, ) cos (2 nfct  +  <pt)  +  n(t) 

i= 1  y=0 

where  as  in  [1]: 

i  =  the  adjacent  cells  /  =  1,2. . .,6 

ij  =  mobile  user  or  channel  j  in  adjacent  cell  i 

Kt  =  the  number  of  active  channels  in  adjacent  cells  i 

Rt  =  Nakagami  fading  random  variable  for  signals  from  adjacent  cells 

Ptj  =  Lognonnal  random  variable  representing  the  average  power  received  from 
the y-th  channel  in  adjacent  cell  i 

by(t)  =  the  information  signal  for  the  y-th  user  channel  in  adjacent  cell  i 
WM)  =  Walsh  function  for  the  y-th  user  channel  in  adjacent  cell  i 
ct{t )  =  PN  spreading  signal  for  the  adjacent  cell  i 
fc  =  the  carrier  frequency  of  the  signal 

tj  =  the  time  delay  from  adjacent  cell  i,  relative  to  the  time  delay  from  the  center 
cell  base  station 
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(p.  =  the  phase  delay  from  adjacent  cell  i,  relative  to  the  phase  delay  from  the 


center  cell  base  station 

The  power  R  can  be  thus  defined  by  consulting  Equation  (2. 12)  as  follows: 


,  f_A 

9  LH(Dt)Xt 


(2.14) 


where 


fv  =  the  factor  used  to  adjust  the  power  in  the y'-th  channel  in  adjacent  cell  i 
Lh  ( D. )  =  the  median  path  loss  using  the  Hata  Model  at  a  separation  distance  Dj 
Z)  =  the  distance  separating  the  receiver  from  the  base  station  in  adjacent  cell  i 
Xi  =  Lognonnal  random  variable  /\  ( 0.  (JiUi ) 

Using  the  received  signal  defined  by  Equation  (2.13),  the  performance  of  the  DS- 
CDMA  system  in  a  Nakagami-fading  channel  with  Lognonnal  shadowing  can  be  defined. 


D.  SUMMARY 


In  this  chapter,  by  taking  into  account  a  number  of  concepts  from  [1],  a 
Nakagami-Lognormal  channel  model  was  developed  which  incorporates  both  large-scale 
and  small-scale  propagation  effects  into  a  single  model.  The  extended  Hata  Model 
predicts  the  median  path  loss  used.  Applying  this  type  of  channel  model  to  the  signal, 
which  is  transmitted  from  the  base  station,  we  end  up  in  a  formula  found  in  Equation 
(2.13)  that  gives  the  received  signal  to  a  mobile  user’s  receiver.  This  received  signal  will 
be  used  to  analyze  the  performance  of  the  DS-CDMA  in  the  following  chapters. 


18 


III.  DS-CDMA  PERFORMANCE  IN  NAKAGAMI-M- 
LOGNORMAL  FADING  CHANNEL 

In  Chapter  II,  a  N a k aga m  i  -m -  Lo g n o rm a  1  channel  model  was  defined,  and  using 
the  same  procedure  as  in  [1],  it  resulted  in  the  received  signal  by  a  mobile  user  on  the 
forward  channel  of  a  DS-CDMA  cellular  system  in  a  Nakagami-m-Lognormal  fading 
environment.  In  this  chapter,  by  using  Forward  Error  Correction  (FEC)  Coding  to  the 
forward  signal  (using  convolutional  codes),  the  coded  probability  of  error  is  computed 
and  the  variation  with  different  values  of  the  Nakagami-m  parameter  can  be  seen.  As 
previously  mentioned,  Nakagami-m  distribution  is  a  more  general  distribution  for 
describing  fading  environments  depending  on  the  value  of  the  m  parameter,  which 
describes  the  severity  of  the  fading  environment. 

The  analysis  which  occurred  in  [1]  considered  that  mobile  user  one,  who  is  the 
receiving  mobile  user,  is  located  in  the  least  optimum  position  within  the  cell,  which  can 
be  any  corner  of  the  hexagon  as  depicted  in  Figure  3.1. 


Figure  3.1.  Mobile  User  One  in  the  Seven-Cell  Cluster. 
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The  distances  from  adjacent  cell  base  stations,  D. . ,  to  mobile  user  one  for  the 
worst-case  distance  d  from  the  base  station  in  the  center  cell  are  given  in  [1]  as  follows: 


D 


=  < 


d, 

2d, 


i  =  4,5 
i  =  3,6 


yfld,  1  =  1,2 


where  d  is  the  length  of  each  side  of  the  hexagon  implying  that  the  hexagon  is  subdivided 
into  six  isosceles  triangles.  Distances  D.  are  used  in  order  to  compute  the  Hata  median 

path  loss  for  the  transmitted  signals  from  base  stations  of  the  i  cell  which  was  previously 
mentioned  as  responsible  for  the  co-channel  interference  in  the  received  signal.  Distance 
d  is  used  for  computing  the  median  path  loss  of  the  information  signal  transmitted  from 
the  center  cell’s  base  station. 


A.  CODED  BIT  ERROR  PROBABILITY 


In  order  to  decode  our  information  bit  stream  the  Viterbi  algorithm  with  soft 
decision  decoding  was  used.  It  is  also  assumed  that  the  all  zero  sequence  (bx(t)  =  1  for  all 
t)  was  transmitted.  The  first  event  error  probability  will  be  used  in  order  to  compute  an 
upper  bound  in  the  coded  bit  error  probability.  The  first  event  error  probability  is  the 
probability  where  at  a  node  B  in  the  trellis  diagram  one  of  the  paths  that  merges  with  the 
all-zero  path  (correct  path)  for  the  first  time  and  the  merging  path  has  a  metric  that 
exceeds  the  all-zero  path  metric  [1]  and  [6].  When  this  situation  occurs,  the  decoder 
discards  the  correct  path  and  bit  errors  ensue.  This  probability  of  error  in  the  pair  wise 
comparison  of  these  two  paths  that  differ  in  d  bits  is  from  [1],  and  [6]  given  by: 


(3.1) 
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where,  Q  denotes  the  Q  function  and  zd  is  the  value  assumed  by  the  new  random 
variable  Zd,  which  is  the  sum  of  d  Nakagami-m  square -Lognormal  random  variables 
defined  as  follows: 

Zj  =  RjXi  (3.2) 

/  =  1 

where 

2 

R,  =  the  Nakagami  square  random  variable, 
d  =  distance  from  the  all  zero  path  in  the  Trellis  diagram, 

Xi=—  is  another  Lognonnal  random  variable  Xi~  K(Q,AcrdB\  since  X,  is 

xi 

our  initial  Lognonnal  random  variable,  X ,  ~  A  ( 0,  Act dB ) .  The  index  /  in  (3.2)  runs  over 
the  set  of  d  bits  in  which  the  two  paths  differ. 


Additionally  we  have  from  [1]: 


exp 


A2  a. 


dB 


a  = 


3  N 


j_  y1  fa 

M  M  f  LH(Dt)  2 Ec 


'y  ^H(d)_+  N0 


(3.3) 


It  has  been  assumed  that  Ec  =  (k0  / n0)Eh  is  the  coded  bit  energy, Eb  =  fxP,T / LH  (d)  is  a 
baseline  received  bit  energy  with  no  fading  or  shadowing  and  n0,  k0  are  the 
characteristics  of  the  encoder  that  is  being  used.  The  use  of  an  (k0,n0)  encoder  means 
for  every  k0  information  bits,  the  encoder  produces  n0  coded  bits  and  that  the  code  rate 
k 

is  R  =  — .  In  order  to  keep  the  same  bit  rate  to  the  system,  the  bit  duration  must  be 

n0 

reduced  as  T  =  TR  =  T  ^"/  [1], 

/  n0 
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It  can  easily  be  said  that  if  the  pdf,  pz  (zd )  is  known,  the  conditioning  of  the 
first-event  error  probability  on  the  new  random  variable  Zd  can  be  removed,  by 
integrating  across  its  pdf  as  follows: 


P2(d)  ^  P2 (d)  zdPzd{zMzd 


f  i — A 
z 


=11  Q  pz  PzAz*)dz, 


(3.4) 


V  ’  ~  J 


However,  in  a  particular  convolutional  code  there  could  be  a  large  number  of  paths  with 
different  distances,  which  could  compete,  with  the  all-zero  path  at  node  B.  That  is  why  it 
is  possible  to  compute  an  upper  bound  on  the  first-event  error  probability.  One  way  to 
calculate  an  upper  bound  on  the  probability  of  bit  error  Pe  is  by  using  the  total  number  of 

information  bit  errors,  f3d,  which  is  determined  by  selecting  a  path  of  distance  d  from  the 
all-zero  path  [6] 

P,di  AA(d)  (3.5) 

^  d=d,  „ . 


Thus,  the  end  result  is  a  formula  which  makes  it  possible  to  calculate  the  probability  of 
error,  Pe ,  for  the  coded  DS-CDMA  cellular  system  in  the  N ak aga m i -m - Lo g n o rm a  1 

channel  depending  on  the  distribution  of  the  new  random  variable  Zd  .  The  calculation  of 
this  pdf  takes  place  in  the  next  Section. 


B.  APPROXIMATING  THE  SUM  OF  MULTIPLICATIVE  NAKAGAMI- 
SQUARE-LOGNORMAL  RANDOM  VARIABLES 


Since  the  new  random  variable  Zd  is  defined  by  Equation  (3.2),  it  would  be 
possible  to  compute  its  probability  density  function  if  the  pdf  for  one  single  variable 
W  =  R2 X  was  known,  which  represents  a  Nakagami-zu-square -Lognormal  random 
variable.  Since  Zd  is  the  sum  of  d  Nakagami-m-square-Lognormal  variables,  it  is  thus 
possible  to  predict  its  pdf  by  convolving  pw{w)  with  itself  d  times.  Additionally,  in 
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order  to  find  the  upper  bound  on  the  probability  of  bit  error  as  predicted  by  Equation 
(3.5),  a  series  of  P2(d)  for  d  =  dfree,  dfrec+ 1,  dfree  +  2,...  must  exist.  Using  the  first  five 

tenns  we  get: 

1  dtiee+4 

p,<-tZ  AW) 


In  our  analysis,  which  is  similar  to  that  occurring  in  [1],  convolutional  encoders  will  be 
used  with  Rcc  =  and  constraint  lengths  form  7  to  9,  which  usually  have  a  dfree  ranging 


from  10  to  12. 


The  sum  of  Lognormal  random  variables  can  be  approximated  as  another 
Lognormal  random  variable  upon  setting  up  a  model  with  suitable  parameters  [1].  There 
are  a  lot  of  different  ways  in  order  to  find  the  Lognormal  model  parameters.  One  of  them 
is  the  cumulant  matching  approach  in  which  the  model  parameters  are  being  computed  to 
make  the  Cumulative  Distribution  Function  (CDF)  of  the  model  Lognormal  random 
variable  fit  the  CDF  of  the  sum  of  Lognormal  random  variables  (as  simulated  using 
Monte  Carlo  methods).  After  defining  the  model  parameters  by  using  one  of  the  methods 
of  [7],  the  approximated  distribution  is  known.  The  other  pdf  to  determine  is  the  pdf  of 
the  sum  of  d  Nakagami-m  square  random  variables,  which  will  be  examined  more  closely 
later  in  this  Section.  Furthermore,  knowing  the  pdf  of  the  sum  of  Lognonnal  variables 
(approximated  as  another  Lognormal  random  variable)  and  the  pdf  of  the  sum  of  d 
Nakagami-zzz  square  random  variables,  the  sum  of  d  Nakagami-nz-square-Lognormal 
random  variables  will  be  approximated  with  Zd  as  a  multiplicative  Nakagami-zzz-square- 
Lognormal  random  variable. 


1.  Defining  the  Model,  Zd 


It  is  known  that  the  Nakagami-zzz  random  variable  has  a  probability  density 
function  that  is  given  by: 
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2  ( m 


fR,{r,)=  r(m)tQ 


exp  — ,r,>  0 


elsewhere 


where  T(m)  is  the  gamma  function,  defined  by  the  integral 


F ( in )  =  [  xm  1  exp {-x)dx,m  >  0 

J  0 


£2  1 
Varfi?,2)  “  2 


q  =  E{f?r  j 


In  Appendix  III-A  it  is  shown  that  the  probability  density  function  and  the 
characteristic  function  of  a  Nakagami-m  square  random  variable  are  given  by: 


/  ,  \  =  2(„-i)  exp  -- 

Jk‘{!  >  r (m)  1  { 


'F  2(j«)  = 


and  that  the  probability  density  function  of  the  sum  of  d  Nakagami-m  square  random 
variables  is  given  by: 


/  \  md 

f  (y  )  =  M  ZM 
Jy-^yid’  [qJ  T(t. 


md  md- 1  x 

y2d  |  w2d 

r  {md)  ^  Q 


The  model  Z.  can  now  be  defined  in  terms  of  its  two  factors,  which  are  defined  as 


follows: 


^  =  y„xm 
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where  Yld  is  a  random  variable  whose  pdf  is  given  by  Equation  (3.9)  and  it  is  the  sum  of  d 

d 

Nakagami-/?z  square  random  variables,  Y2d  =  ^  R~  and  Xz  is  the  lognormal  component 

1=1 

of  Zd,X.  ~  A(jUz,crz) .  The  parameters  jiz  and  a.  are  yet  to  be  specified. 


The  quantity  E  jf?,2  j  =  £2  =  1  has  been  normalized.  For  the  moments  of  Yld : 


E  {Yld]  =  dE{Rj}  =  dQ.  =  d 

E \Y2d2}  =  Qrd2  +  £22  —  =  d2  +  — 
m  m 

Var{FM}  =  Q!-  =  - 
m  m 


(3.10) 


Assuming  that  Yld  and  X.  are  independent  random  variables,  the  expectation  for  the 
model  is: 

E{z„}  =  E{F„X)  =  E{F„}E{X)  (3.11) 

From  [1],  it  is  known  that 

(  fT2'] 

E{X,}  =  exp  (3-12) 


Substituting  the  moments  of  the  two  independent  random  variables  Yld  and  X_  in 
Equation  (3. 1 1),  we  get  for  the  expectation  of  the  model: 

EjZrf}  =  Yld exp(//_  +07  / 2)  (3.13) 

The  second  moment  of  Zd  can  be  defined  from  the  second  moments  of  Y2d  and 
X.  as  follows: 

E{z;}  =  E{r2X!}  =  E{ri}E{^;}  (3.14) 

From  [1], 

E{x?}  =  exp(2/tr+2f7;)  (3.15) 
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Substituting  the  values  of  E j Y^d j  and  e{x2}  from  Equations  (3.10)  and  (3.15) 
respectively,  we  get  for  E  jzj  j  : 


E{Z;}  =  Q3 


d2  + 


d_ 
m  J 


exp(2//_+2cr2) 


(3.16) 


knowing  the  first  and  second  moment  of  Zd  the  variance  of  Zd  can  be  computed  as 
follows: 


Var {z,}-- 


( 


Q2 


d 


d2  + 

y  mj 


exp  (2 juz  +  2 cr; )  -  d2  exp(2//_  +  (J  2)Q 


(3.17) 


However,  the  parameters  flz  and  a:  for  the  random  variable  Xz  have  not  yet 
been  determined  in  order  to  approximate  Zd  as  a  Nakagami-nz-square  Lognormal 
random  variable.  In  order  to  calculate  these  two  parameters  we  relate  the  mean  and  the 
variance  of  the  model  Zd  defined  by  Equations  (3.13)  and  (3.17)  respectively  with  a 
scaled  version  of  mean  and  variance  of  the  single  Nakagami-/«-square  Lognormal 
random  variable  R2 X . 


The  previous  analysis  is  similar  to  that  taking  place  in  [1]  in  which  the  random 
variable  R  is  a  Rayleigh  random  variable,  special  case  of  the  Nakagami-m  random 
variable  (m  =  1).  As  in  [1],  in  simulations  it  can  be  seen  that  the  mean  of  Zd  is 
proportional  to  d  times  the  mean  of  the  original  single  Nakagami-zu-square-Lognormal 
random  variable.  Based  on  that,  the  mean  of  Zd  can  be  equated  with  d  times  the  mean  of 

R2X  scaled  by  a  factor  of  fx  as  follows: 


E{zrf}  =  fdE{R2X }  =  /</E{R2}E{X}  =  fdQex  p 


(  CT2^ 

(  2y 

/T 

=  fd  exp 

°x 

l  2  J 

l  2  ) 

(3.18) 


where  a2  =  dd a] H  is  the  parameter  from  the  original  Lognormal  shadowing  random 
variable  X  ~  A  ( 0,  cr2- ) .  Using  the  Equations  (3.13)  and  (3.18): 
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(3-19) 


The  analytical  computation  for  jUz  is  given  in  Appendix  III-C.  Observing  that  the 
variance  of  Zd  is  proportional  to  d  times  the  variance  of  the  Nagakami-/?z-square- 
Lognormal  random  variable;  the  variance  of  the  model  can  be  equated  with  d  times  the 
variance  of  R2X  scaled  by  a  factor  of  f2  as  follows: 

Var{zJ  =  /2JVar{i?2X}  (3.20) 

The  variance  of  R2 X  is  calculated  in  Appendix  III-B  and  it  is  given  as  follows: 

VarjiCx}  =  expfrr2)  Q2 ea' - Q2  (3.21) 

1  1  V  L  m 

By  using  Equations  (3.20)  and  (3.21),  <72  can  be  solved  as  follows: 

Gz  =  In  —  +  —  en’:  -  +  d  -  In  ( 1  +  md )  +  ln(/;z)  (3 .22) 

fm  f 

Analytical  computation  for  a\  is  given  in  Appendix  III-C.  The  scaling  factors  f\ 
and  f2  vary  with  the  values  of  adB  and  the  Nakagami-zu  parameter,  and  have  been 
selected  in  the  way  that  the  CDF  of  the  model  Zd  best  fits  the  CDF  of  Zd  for  d  = 
8, 9...,  16.  The  value  of  odB  =  7  will  be  used  to  examine  the  perfonnance  of  the  model 
Zd  for  the  following  values  of  the  Nakagami-nz  parameter:  m  =  0.5,  0.75,  1,  1.5,  2.  Each 
value  of  in  corresponds  to  a  different  fading  environment.  Table  3.1  shows  the  values  of 
f\  and  f2  for  each  value  of  odB  as  determined  by  simulation,  for  m  =  0.5. 
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The  model  parameters  have  been  specified,  so  it  is  now  possible  to  define  the  pdf 
of  the  model  Zd  using  the  marginal  pdfs  of  its  components,  which  are  Y2d  and  X _ .  The 


quantity  Zd  —  YldX  _  has  been  defined  and  X,  —  X  can  be  fixed  so  that  Zd 
The  conditional  pdf  of  Zd  I  v  is  given  as  follows: 


lx 


However,  the  pdf  pY  /  of  Yld  is  given  by  Equation  (3.9)  so  (3.23)  becomes: 


f  z\md~x 


%(ZW=7 


/  \  md 

1  m  ' 
\£lj 


T{md) 


exp 


^  m  z^ 

Q  Xy 


1 

md  i  md- 1 

1  Z 
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m  z' 

X 

T(md)  xmd~l  P 

V 

Q  X  y 

/  \  md 

1  m  ' 

T{md)xmd 


md— 1  „  „ 

■z  exp 


m  z 


Q.  xj 


->0 


The  conditional  of  Zd  on  X,  —  x  can  be  removed  as  follows: 


=  xy2d- 


(3.23) 


(3.24) 
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PiSz)=l 


Pn,  P<z(x)dx 


=r 

Jo 


m 

£2 


v  md 


T(md)x 


md 


md—\ 
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m  z 

£2  x 


xx/2 Ha, 


x  exp 
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md 


m 

£2  j 
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Jo  ymd+l 


exp 


i 

-zm  (in  x~nz) 


x£2 


2cr 


(3.25) 


dx ,  z  >  0 


2 

where  jU,  and  <7~z  have  been  defined  by  Equations  (3.19)  and  (3.22)  respectively. 

A  distribution  for  the  Nakagami-zn-square-Lognormal  random  variable  has  now 
been  defined,  which  approximates  the  distribution  of  a  sum  of  d  Nakagami-zu-square- 
Lognormal  random  variables.  In  the  next  Section,  some  examples  of  the  model  will  be 
examined  and  the  perfonnance  results  using  the  model  with  the  results  predicted  by  the 
Monte  Carlo  simulations  will  be  compared. 


2.  Testing  the  Model  Zd 

In  this  Section,  the  model’s  applicability  to  bit  error  analysis  of  DS-CDMA  with 
forward  error  correction  coding  will  be  demonstrated.  Using  the  first  five  terms  of  the 
union  bound  from  Equation  (3.5)  we  will  approximate  the  probability  of  bit  error.  The 
value  of  lognormal  shadowing  will  be  taken  to  be  adB  =  7  and  a  convolutional  encoder 

with  a  rate  of  and  a  constraint  length  v  =  8  will  be  used.  The  function  P2(d)  will  be 

approximated  using  the  model  for  d  =  10  through  14  which  are  critical  values  of  d  for  this 
type  of  encoder  with  constraints  lengths  7  and  8.  MATLAB  is  the  programming  software 
used  to  obtain  my  results.  First  of  all,  the  model  will  be  checked  and  compared  with  the 

one  developed  in  [1]  by  giving  the  Nakagami-/«  parameter  the  value  of  m  =  1  (which  is 
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the  Rayleigh  fading  case),  if  it  reproduces  the  same  results.  Figures  3.2  through  3.6  show 
histograms  for  Zd  for  udB  =  7  and  Nakagami-m  parameter  m  =  1  with  an  over  plot  of  the 

pdf  of  our  model.  The  histograms  for  Zd  were  developed  by  generating  d  sets  of 

Nakagami-square-Lognormal  distributed  data  consisting  of  400,000  independent  samples 
per  set  and  summing  over  the  d  sets  similar  to  those  in  [1].  In  these  Figures  concerning 
the  histograms  part,  the  x  axis  represents  the  possible  values  of  the  random  variable  we 
are  trying  to  model  and  the  y  axis  is  the  number  of  occurrences  of  the  possible  values  of  x 
within  each  bin.  For  the  pdf  part  the  v  axis  is  the  value  of  the  pdf  normalized  by  the 
number  of  samples  we  generate  and  the  binwidth. 
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Figure  3.2.  Histogram  of  Z10  and  the  PDF  for  Z10  for  adB  =  7 ,  d  =  10  and  m  =  1. 
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Figure  3.3.  Histogram  of  Zn  and  the  PDF  for  Zn  for  odB  =7  ,d  =  11  and  m  =  1. 


35 


x  104 


0  20  40  60  80  100  120  140  160  180  200 


Figure  3.4.  Histogram  of  Z12  and  the  PDF  for  Z12  for  udB  =  7  ,  d  =  12  and  m  =  1. 
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Figure  3.5.  Histogram  of  Z13  and  the  PDF  for  Zu  for  adB  -1  ,d  =  13  and  m  =  I . 
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Figure  3.6.  Histogram  of  Z14  and  the  PDF  for  Z14  for  udB  =7  ,d=l4  and  m  =  1. 
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Comparing  these  results  with  the  corresponding  results,  which  were  obtained  in 
[1],  it  can  be  seen  that  the  model  Zd  has  almost  the  same  results  as  the  model  in  [1]  for 
the  Rayleigh  case. 

Figures  3.7  through  3.11  show  the  histograms  for  Zd  for  a dB  =  7  and  the  pdf  of 
the  model  Zd  for  m  =  0.5  and  for  d  =  10  through  <7=14. 
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Figure  3.7.  Histogram  of  Z10  and  the  PDF  for  Z10  for  adB  =7  ,d  =  10  and  in  =  0.5. 
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Figure  3.9.  Flistogram  of  Z12  and  the  PDF  for  Z12  for  udB  =7  ,d  =  12  and  m  =  0.5. 
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Figure  3.10.  Histogram  of  Z13  and  the  PDF  for  Z13  for  <JdB  =7  ,d  =  13  and  m  =  0.5. 
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Figure  3.11.  Histogram  of  Z14  and  the  PDF  for  Z14  for  adB  =7  ,d  =  14  and  m  =  0.5. 
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Figures  3.12  through  3.16  show  histograms  for  Zd  for  crdB  -  7  and  m  =-  0.75  and 


an  over  plot  of  the  pdf  of  our  model  for  d  =  10  through  <7=14. 
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Figure  3.12.  Histogram  of  Z10  and  the  PDF  for  Zi0  for  <JdB  =7  ,d  =  10  and  m  =  0.75. 
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Figure  3.13.  Histogram  of  Zn  and  the  PDF  for  Zn  for  a dB  =7  ,d  =  11  and  m  =  0.75. 
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Figure  3.14.  Histogram  of  Z12  and  the  PDF  for  Zl2  for  adB  =7  ,d  =  12  and  m  =  0.75. 
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Figure  3.16.  Histogram  of  Z14  and  the  PDF  for  Z14  for  adB  =7  ,d  =  14  and  m  =  0.75. 

Figures  3.17  through  3.21  show  histograms  for  Zd  for  odB  =7  and  m  =  1.5  and 
an  over  plot  of  the  pdf  of  our  model  for  d  =  10  through  d=  14. 
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Figure  3.17.  Histogram  of  ZI0  and  the  PDF  for  Z10  for  odB  =7  ,d  =  10  and  m=  1.5. 
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Figure  3.19.  Flistogram  of  Zn  and  the  PDF  for  Z12  for  udB  =7  ,d  =  12  and  m  =1.5. 
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Figure  3.20.  Histogram  of  Z13  and  the  PDF  for  ZI3  for  adB  =7  ,d  =  13  and  m  =1.5. 
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Figure  3.21.  Histogram  of  Z14  and  the  PDF  for  Z14  for  adB  =7  ,d  =  14  and  m  =  1.5. 

Figures  3.22  through  3.26  show  histograms  for  Zd  for  adB  =  7  and  m  =  2  and  an 
over  plot  of  the  pdf  of  our  model  for  d  =  10  through  d  =  14. 
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Figure  3.22.  Histogram  of  Z10  and  the  PDF  for  Z10  for  adB  =  7  ,  d  =  10  and  m  =  2. 
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Histogram  of  Zn  and  the  PDF  for  Zn  for  a dB  =7  ,d  =  11  and  m  =  2. 
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Figure  3.23. 


Figure  3.25.  Histogram  of  Z13  and  the  PDF  for  Z13  for  crdB  =7  ,d  =  13  and  m  =  2. 


Figure  3.26.  Histogram  of  Z14  and  the  PDF  for  Z14  for  odB  -  7  ,  d  =  14  and  m  =  2. 


43 


The  next  question  is  how  well  can  the  model  Zd  approximate  Zd  in  calculating 
the  first-event  error  probability,  which  is  given  by  (3.4).  This  answer  will  be  given  by 
comparing  the  result  predicted  by  the  model  Zd  to  that  one  which  simulates  the  integral 
of  (3.4)  by  generating  d  independent  samples  form  the  N akaga i m -in - sq uare- Logn o rm a  1 
distribution  and  summing  them  in  order  to  form  one  realization  for  Zd  and  one 

realization,  px,  for  P:(d) .  This  process  is  repeated  95,000  times  and  gives  our 
estimation,  p.,  for  P2(d)  as  follows: 


P  = 


1 

95,000 


95,000 


I  Pi 


(3.26) 


First  of  all,  by  using  the  following  values,  an  attempt  will  be  made  to  replicate  the 
results  obtained  in  [1]  using  the  model  Zd  as  follows:  in  =  1  (Rayleigh  case),  adB  =  7  , 
number  of  sectors  6,  60  users  per  adjacent  cell  and  d  =  10  through  13.  Figures  3.27 
through  3.36  show  how  well  the  model  Zd  is  useful  in  computing  P2(d) and  comparing  it 
to  that  predicted  by  Monte  Carlo  simulations  of  the  integral  of  (3.4). 
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Figure  3.27.  First  Event  Error  Probability  (1 0)  and  P2{\  1)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  1 . 


Figure  3.28.  First  Event  Error  Probability  73,(12)  and  73(13)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m=  1 . 
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Figure  3.29.  First  Event  Error  Probability  /3(H))  and  P2{\  1)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  0.5. 


Figure  3.30.  First  Event  Error  Probability  P2(\ 2)  and  /3 ( 1 3 )  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  0.5. 
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Figure  3.31.  First  Event  Error  Probability  P2  (1 0)  and  P2  (1 1)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  0.75. 


Figure  3.32.  First  Event  Error  Probability  73,(12)  and  73(13)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  0.75. 
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Figure  3.33.  First  Event  Error  Probability  /3(H))  and  P2(  1 1)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  1.5. 


Figure  3.34.  First  Event  Error  Probability  P2 (1 2)  and  /3 ( 1 3 )  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =1.5. 
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Figure  3.35.  First  Event  Error  Probability  /3(H))  and  P2{\  1)  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  2. 


Figure  3.36.  First  Event  Error  Probability  P,(  12)  and  /3 ( 1 3 )  with  60  Users  Per  Cell 

and  60°  Sectoring,  for  m  =  2. 
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Observing  Figures  3.27  through  3.36,  it  can  be  seen  that  the  model  Zd 
approximates  the  simulated  results  well  for  P2(d )  for  the  values  of  in  =  0.5,  0.75,  1.5  and 
2.  Using  these  results  for  P2(d)  and  adding  them  to  (3.5),  it  is  possible  to  compute  the 
modeled  probability  of  bit  error  for  the  first  five  terms  in  the  union  bound.  Figures  3.37 
through  3.41  show  the  comparison  of  the  modeled  and  simulated  probability  of  bit  error 
using  different  values  of  m.  To  obtain  these  results,  a  rate  of  Vi  convolutional  encoder 
was  incorporated  with  a  length  of  v  =  8.  The  average  received  SNR  per  bit  for  the 
Nakagami-Lognormal  channel  is  from  [1]: 


since  E{R2]  =  £2  =  1  has  been  normalized. 


Figure  3.37.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  with  m  =  1 
(Rayleigh  Case)  and  Lognormal  Shadowing  ( adB  =  7)  using  Six  Sectors  and  a  Rate  Vi 

Convolutional  Encoder  with  v  =  8. 
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For  the  case  of  v  =  8,  the  modeled  and  simulated  results  of  P2(d)  were  obtained 
for  d  =10  through  14  since  dfree  =  10  and  j3l0  =  2,  j3n  =  22,  (3n  =  60,  j3n  =  148,  (3U  =  340 

[1].  By  selecting  the  Nakagami  parameter  to  be  m  =  1,  the  channel  becomes  a  Rayleigh- 
Lognormal  channel. 

Comparing  the  modeled  and  simulated  probability  of  bit  error,  it  can  be  observed 
that  the  modeled  and  simulated  results  do  not  differ  much.  More  specifically,  at  the 
range  10  to  15  dB  (probability  of  bit  error  10°  to  10”4)  the  differences  between  the 
modeled  and  simulated  results  are  quite  small.  This  particular  range  of  SNR  is  a  practical 
one  in  which  DS-CDMA  cellular  systems  operate.  Differences  between  the  modeled  and 
simulated  results  become  notable  when  the  interference  floor  is  in  a  range  between  10'4 
and  10'7.  The  difference  in  the  probability  of  bit  error  as  the  number  of  users  per  cell 
increases  from  10  to  60  should  also  be  mentioned.  Channel  bit  error  probability  becomes 
worse  as  the  number  of  users  in  a  cell  increases  due  to  co-channel  interference. 

Figure  3.38  is  similar  to  Figure  3.37,  that  compares  the  modeled  and  simulated 
probability  of  bit  error  but  for  a  different  value  of  m.  In  this  case,  the  Nakagami 
parameter  m  equals  to  0.5.  This  type  of  fading  is  more  severe  than  Rayleigh  fading.  The 
differences  in  the  fading  environment  as  the  value  of  in  varies  will  be  examined  later  in 
this  chapter.  As  for  the  differences  in  the  modeled  and  simulated  results,  the  same 
comments  from  the  previous  figure  apply,  but  in  this  case,  the  probability  of  bit  error  is 
worse  when  comparing  to  the  previous  one  for  all  the  range  values  of  SNR. 
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Probability  of  bit  error 


Figure  3.38.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  with  m  = 
0.5  and  Lognormal  Shadowing  (<7dfj  =  7)  using  Six  Sectors  and  a  Rate  V2  Convolutional 

Encoder  with  v  =  8. 
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Figure  3.39.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  with  m  = 
0.75  and  Lognormal  Shadowing  ( udB  =  7)  using  Six  Sectors  and  a  Rate  Vi  Convolutional 

Encoder  with  v  =  8. 


In  the  case  of  the  Nakagami  parameter  m  =  0.75,  the  modeled  and  simulated 
results  for  the  bit  error  probability  are  very  similar  in  the  range  of  SNR  from  10  to  15  dB. 
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Figure  3.40.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  with  m  = 
1.5  and  Lognormal  Shadowing  (<JdH  =  7)  using  Six  Sectors  and  a  Rate  lA  Convolutional 

Encoder  with  v  =  8. 


In  Figure  3.40,  it  can  be  seen  that  the  differences  between  the  simulated  and 
modeled  probability  of  error  are  noteworthy  and  greater  than  the  previous  differences  in 
the  SNR  range  of  20  to  30  dB  but  they  correspond  to  the  interference  floor  from  1 0  4  to 
10'7  of  bit  error  probability. 
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Figure  3 .4 1 .  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  with  m  =  2 
and  Lognormal  Shadowing  ( crdB  =  7)  using  Six  Sectors  and  a  Rate  Vi  Convolutional 

Encoder  with  v  =  8. 


Figure  3.41  depicts  that  our  modeled  and  simulated  results  fit  each  other  in  the 
SNR  range  10  to  15  dB  in  the  case  of  the  Nakagami  parameter  m  =  2  which  corresponds 
to  Ricean  fading.  There  is  LOS  path  between  the  base  station  and  the  mobile  user. 
Appendix  III-D  shows  similar  results  for  adB  -  4  .  The  objective  of  the  next  Section  is  to 

use  these  results  to  explore  the  performance  analysis  of  our  system  in  a  Nakagami  Fading 
and  Lognormal  Shadowing  environment. 


C.  BIT  ERROR  ANALYSIS  OF  DS-CDMA  WITH  FEC 


In  the  previous  Section,  a  model  was  created  which  is  a  useful  tool  for 
approximating  the  distribution  for  a  sum  of  multiplicative  Nakagami  square-Lognormal 
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random  variables.  In  this  Section,  by  using  this  approximation,  the  performance  of  DS- 
CDMA  with  FEC  will  be  analyzed  for  the  different  values  of  the  Nakagami  parameter  m. 
Figures  3.42  through  3.47  show  these  results  for  20,  40,  and  60  users  per  adjacent  cell 
respectively  for  both  modeled  and  Monte  Carlo  simulated  results.  In  each  case,  the 
Nakagami  parameter  m  varies  for  the  specific  values  of  m  =  0.5,  in  =  0.75,  in  =  1,  m  =  1.5 
and  m  =  2. 


SNR  per  bit  (dB) 


Figure  3.42.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  with  FEC  in 
Various  Fading  Conditions  with  20  Users  Per  Cell,  Employing  a  60°  Sectoring  ( Rcc=  1/2 
and  v  =  8)  and  Lognormal  Shadowing  of  adB  =  7  . 
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Probability  of  bit  error 


SNR  per  bit  (dB) 

Figure  3.43.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  with  FEC  in  Various  Fading  Conditions  with  20  Users  Per  Cell,  Employing  a 
60°  Sectoring  (Rcc=  1/2  and  v  =  8)  and  Lognormal  Shadowing  of  adB  =  7  . 
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Probability  of  bit  error 


SNR  per  bit  (dB) 

Figure  3.44.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  with  FEC  in 
Various  Fading  Conditions  with  40  Users  Per  Cell,  Employing  a  60°  Sectoring  (Rcc  =  1/2 
and  v  =  8)  and  Lognormal  Shadowing  of  adB  -  7  . 


58 


Probability  of  bit  error 


SNR  per  bit  (dB) 


Figure  3.45.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  with  FEC  in  Various  Fading  Conditions  with  40  Users  Per  Cell,  Employing  a 
60°  Sectoring  (Rcc=  1/2  and  v  =  8)  and  Lognormal  Shadowing  of  adB  =  7  . 
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Probability  of  bit  error 


SNR  per  bit  (dB) 

Figure  3.46.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  with  FEC  in 
Various  Fading  Conditions  with  60  Users  Per  Cell,  Employing  a  60°  Sectoring  (Rcc  =  1/2 
and  v  =  8)  and  Lognormal  Shadowing  of  adB  -  7  . 
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Figure  3.47.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  with  FEC  in  Various  Fading  Conditions  with  60  Users  Per  Cell,  Employing  a 
60°  Sectoring  (Rcc=  1/2  and  v  =  8)  and  Lognormal  Shadowing  of  <7dB  =  7  . 

From  Figures  3.42,  through  3.47,  it  can  be  seen  that  when  our  channel  operates  in 
a  fading  environment  where  m  =  0.5,  it  performs  worse  than  it  does  for  the  other  four 
values  of  m.  On  the  other  hand,  when  m  =  2,  the  case  of  Ricean  Fading,  a  LOS  path 
occurs  and  our  channel  performs  better  than  it  does  for  the  other  types  of  fading.  An 
average  fading  condition  between  them  is  the  Rayleigh  Fading  environment  where  m  =  1. 
There  is  no  great  difference  in  the  perfonnance  of  our  system  when  the  value  of  m  goes 
from  1.5  to  2.  It  should  also  be  noted  how  the  interference  floor  varies  as  the  number  of 
users  increases  per  adjacent  cell.  In  the  case  of  Figure  3.42  with  20  users  per  adjacent 
cell,  the  interference  floor  corresponds  to  a  probability  of  bit  error  from  1 0  4  to  10  \ 

■j 

Furthermore,  for  SNR  =15  dB,  our  system  achieves  a  probability  of  bit  error  from  10'  to 
10'4  for  the  five  types  of  fading  environment.  As  the  number  of  users  per  adjacent  cell 
increases  from  20  to  40,  as  shown  in  Figure  3.44,  the  intercell  interference  increases.  The 
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new  level  now  is  above  10‘4  of  bit  error  probability.  For  the  same  value  of  SNR  (15  dB) 
our  system  can  now  achieve  a  bit  error  probability  of  10'  to  10'  for  the  different  types  of 
fading.  Comparing  the  results  from  Figures  3.42  and  3.44,  our  system  performs  better  in 
Figure  3.42  because  of  the  lesser  amount  of  intercell  interference.  Additionally,  results 
become  worse  if  the  number  of  users  per  adjacent  cell  increases  to  60  as  seen  in  Figure 
3.46.  For  SNR  =15  dB,  a  bit  error  probability  of  10'  to  10'“  can  be  achieved  and  the 
interference  floor  is  above  10‘3  to  10‘4.  Similar  figures  are  obtained  in  Appendix  III-E  for 
adH  =  4 .  If  the  60°  sectoring  directional  antennas  with  constant  gain  inside  the  60° 

antenna  pattern  and  zero  dB  gain  outside  of  it,  are  replaced  with  120°  antennas  or  with 
one  antenna  (no  sectoring  case),  the  amount  of  intercell  interference  increases  by  a  factor 
of  3  and  6  respectively.  Figures  3.48  through  3.57  show  the  results  obtained  using  the 
model  and  the  Monte  Carlo  simulations  with  six  sectors,  three  sectors  and  one  sector. 
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Figure  3.48.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  Using 
Sectoring  for  m  =  0.5  and  adB  =  7  (Rcc  =  1/2  and  v  =  8). 
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Figure  3.49.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  Using  Sectoring  for  m  =  0.5  and  adB  =  7  (Rcc  =  1/2  and  v  =  8). 
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Figure  3.50.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  Using 
Sectoring  for  m  =  0.75  and  adB  =  7  (Rcc  =  1/2  and  v  =  8). 
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Figure  3.51.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  Using  Sectoring  for  m  =  0.75  and  <JdB  =  7  (Rcc  =  1/2  and  v  =  8). 
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Figure  3.52.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  Using 
Sectoring  for  m  =  1  and  adB  =  7  (Rcc=  1/2  and  v  =  8). 
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Figure  3.53.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  Using  Sectoring  for  m  =  1  and  adB  =  7  ( Rcc  =  1/2  and  v  =  8). 


67 


Probability  of  bit  error 


10 


10  r 


10  r 


10  r 


10  r 


10  r 


10 


0 


20  users  per  cell 


m=1.5 


\N  \ 


‘X 


No  sectoring 


I,  X 

w 


•X  -  -x  -  ■ 


-X-  -  *-  -  X-  -  X- 


\  X 


' 

'  ->  3  sectors 

\ 

X  '*'■  - 

\ 

\ 

\ 

x\  6  sectors 

S 

X 

x 

'x-  - 


x  -  -x-  - 


■X"  '  -X - -K-  -  •  ( 


10  15  20 

SNR  per  bit  (dB) 


25 


30 


35 


Figure  3.54.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  Using 
Sectoring  for  m  =  1.5  and  adB  =  7  (Rcc  =  1/2  and  v  =  8). 


68 


Probability  of  bit  error 


10 


10'1  t 


10'2  r 


10'3  r 


10'4  r 


10’5  r 


10 


V 


&'q  20  users  per  adjacent  cell 

B  SQ 


m=1 .5 


'b  <\ 

"  x  o. 

'cx  No  sectoring 

V©  ^  -O-O'O-- 

\  % 

b\ 

'  © 

b  ''© 

\  \ 

\  ©x 

©  N  'O  3  sectors 


6  sectors 


10 


X>. 


b 


©-. 


©- 


■o--©'0-b 


© 


15  20 

SNR  per  bit  (dB) 


25 


30 


35 


Figure  3.55.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  Using  Sectoring  for  m  =  1.5  and  adB  =  7  (Rcc=  1/2  and  v  =  8). 
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Figure  3.56.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Model  Using 
Sectoring  for  m  =  2  and  adB  =  7  (Rcc  =  1/2  and  v  =  8). 
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Figure  3.57.  Probability  of  Bit  Error  for  DS-CDMA  Predicted  by  Monte  Carlo 
Simulations  Using  Sectoring  for  m  =  2  and  adB  =  7  ( Rcc  =  1/2  and  v  =  8). 


The  degradation  in  performance  is  obvious  when  changing  to  120°  sectoring  or  no 
sectoring  in  all  the  above  cases. 


D.  SUMMARY 


In  this  chapter,  given  the  expression  for  the  upper  bound  on  the  probability  of  bit 
error,  the  performance  of  the  DS-CDMA  cellular  system  operating  in  a  Nakagami-m- 
Lognormal  channel  was  analyzed  under  various  operating  conditions.  However,  in  order 
to  be  able  to  analyze  the  DS-CDMA  cellular  system  performance,  the  sum  of  d 
Nakagami-square-Lognonnal  random  variables  is  approximated  as  a  multiplicative 
Nakagami-square-Lognormal  random  variable  by  the  use  of  the  model  Zd .  Our 
performance  analysis  took  into  account  that  our  intended  mobile  user  was  in  the  worst- 
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case  position  at  the  edge  of  the  center  cell  shown  in  Figure  3.1,  and  that  60°  sectoring 
was  employed  by  using  six  directional  antennas  at  the  base  stations.  Using  this  type  of 
antenna  at  the  base  station,  the  number  of  interferers  decreased  from  6  to  1  in  the  first  tier 
assuming  a  7-cell  reuse  pattern.  Furthermore,  by  varying  the  Nakagami  parameter  m,  it 
was  possible  to  see  how  our  system  performed  in  various  severe  types  of  fading.  The 
worst  case  of  system  performance  was  the  one-sided  Gaussian  fading  environment,  which 
corresponds  to  a  value  of  the  Nakagami-/;/  parameter  of  0.5.  On  the  other  hand,  the  less 
severe  fading  environment  was  Ricean  Fading  (in  =  2)  where  there  is  a  LOS  path  between 
the  transmitting  base  station  and  the  intended  user.  The  improvement  in  our  system’s 
performance  is  obvious  when  changing  from  in  =  0.5  to  in  =  2.  In  the  next  chapter,  more 
realistic  results  will  be  shown  by  assuming  that  the  mobile  user  is  randomly  distributed  in 
the  cell. 
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APPENDIX  III-A.  DEVELOPMENT  OF  THE  PDF  OF  THE 
NAKAGAMI-SQUARE  RANDOM  VARIABLE  AND  THE  PDF  OF 
THE  SUM  OF  INDEPENDENT  NAKAGAMI-SQUARE  RANDOM 

VARIABLES 


It  is  known  that  the  pdf  of  a  Nakagami  random  variable  R,  is  given  as  follows: 


T(m) 


/  \ 

m 

(  2\ 

m  ) 

2m-\ 

—mrt 

— 

rt  exp 

v  Q  J 

(3.28) 


Given  (3.28),  we  can  compute  the  pdf  of  the  Nakagami- square  random  variable  R ,  using 

2  2 

the  transformation  Y  =  R:  .  Then  for  the  pdf  of  R,  we  get: 

1 


fy(y)  ^  ^ j—  Ir,  (  \[y )  J r i  (  Vr  ) 


(3.29) 


and  since  our  random  variable  has  only  positive  values,  Equation  (3.29)  becomes: 


fAy)  =  ^fRl(4y) 


1  (  M  )  I — 2(m-l) 


r(/«)Uf 


& 


exp 


f  -my  ^ 
V  Q  j 


(3.30) 


2 

Substituting  Y with  R~  we  get: 


fAri)=  1 


r(m) 


f  m  3 


2(m— 1) 

rt  exp 


^  -mi)  ^ 
V  Q  J 


(3.31) 


Given  the  pdf  of  the  Nakagami-square  random  variable  we  can  compute  the  characteristic 
function  of  the  Nakagami-square  random  variable  as  follows: 


^y{ j®)=  J  fY{y)™p()COy)dy 

1 


r  ( m ) 


/  oo 

m  r  m_. 

( m  .  3 

—  y  exp 

— -J« 
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w  { 

J 

(3.32) 


dy 
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By  the  use  of  the  following  identity  from  [8]: 


Jjcv  1  exp  (-//*)  bh:  =  Re(/t)  >  0,  Re(v)  >  0 

o  ^ 


(3.33) 


where  r(v)  is  the  gamma  function  defined  by: 


T  ( v)  =  J  xv  1  exp(-,v)bh:,  v  >  0 


(3.34) 


and  make  the  following  substitutions: 


—  -]CO  =  ju  and  v  =  ju  we  get  for  the  characteristic  function  of  the  Nakagami-square 


random  variable  the  following  equation: 


xu  i  -  \  1  f  mT  T(m) 

V  '  r (m  (m  . 


(3.35) 


and  thus: 


(3.36) 


Now  since  R,  are  independent  random  variables  the  characteristic  function  of  the 


u 

random  variable  Yld  =  is  given  as  follows: 
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f  m  3 

_  yQy 


/  \  md 

'  m  . 

- ]&> 

Q  J 


(3.37) 


so  the  pdf  of  the  sum  of  d  Nakagami-square  random  variables  is  given  by: 


fv2d(y2d)- 2n] 


s  md 


exp(-j  coy2d)dco 


m 

Q 


-J  co 


(3.38) 


in 


Substituting  —  =  ax  and  md  =  a2 ,  Equation  (3.38)  becomes: 


1  00  oc  ^ 

fr2d  ( )  =  —  J  - 1 — ^exp(-j coy2d)dco 


(at-jco) 


(3.39) 


Furthermore  we  can  use  the  following  identity  from  [8]: 


,  p>0 


(fi-)x) 


r(v) 


(3.40) 


and  making  the  substitutions  (3  =  ccx,  v  =  a2  and  p  =  y2d ,  we  get  for  fY 


Vld  In  1  r(a2) 


/  md- 1 

/  \  1  J 
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yQ.j 
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F  ( md ) 


exp 


myid 

v  ^  y 


(3.41) 
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APPENDIX  III-B.  VARIANCE  OF  R2X 


The  variance  of  the  product  of  the  Nakagami-square  random  variable  and  the 
Lognormal  random  variable  can  be  computed  as  follows.  Assume  that  A  =  R2  X  ,  then 


E  {A2}  =  E{i?4X2}  =  E{f?4  |e{x2  j 


(3-42) 


Var{f?2x}  =  E{i?4X2}-[E{f?2x} 


(3.43) 


The  characteristic  function  of  R2  is  as  follows: 


VAin) 


f  m  3 


(3.44) 


in 

Q 


-J  O) 


Taking  the  derivative  of  both  sides  in  Equation  (3.44)  with  respect  to  co ,  the  following 
appears: 
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(3.45) 


The  second  derivative  of  the  characteristic  function  is: 
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^;2a  co) 


.  m  /  m  .  ,  .. 

-J  o  (m  +  1)  o_J"  (_J) 


m  —  ( m  +  1 ) 


(3.46) 


It  is  already  assumed  that  e{7?2}  =  Q  =  1  and  it  is  known  that: 


i~ 


m  —  ( m  +  1 ) 

,  J 


(m  +  \)£l2 


(3.47) 


From  [1]  it  is  also  given  that: 


E{zJ}=exp(2<T;) 


(3.48) 


Applying  Equations  (3.47)  and  (3.48)  into  (3.42): 


E{x42  j  =  exp(2<r2) 

1  J  m  v  ' 


(3.49) 


Since  E[i?2x]  =  E|i?"]E{X}  =  Q7  exp(<72  / 2 )  and  Q2  =  1 ,  Equation  (3.43)  becomes: 
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Var{i?2x}  =  ’-^-Q2  exp  (la2 )  -  Q2  exp  (a2 ) 

=  exp(cr)Q2  ^^-exp(cr2)-l  (3.50) 

\_  m 

=  exp  (a2 )  ^-Uxp(cr2)-1 

|_  m 
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APPENDIX  III-C.  MEAN  AND  VARIANCE  OF  THE  MODEL 


Equation  (3.18)  gives  that  E{Zd}  =  fd exp 


V  2  J 


.  Equation  (3.13)  gives  that 


f  2  X 

/x 

{Zd}  =  Qd  exp 

/^z  H - — 

L  2J 

.  Equating  the  last  two  equations  we  get: 
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(3.51) 


Equation  (3.17)  gives: 


Var{z,}  =  £2Jl  rf2  +  ^Jexp(2ft  +2^) 
-d2  exp(2//_  +  2ct2)Q2 


Equation  (3.20)  gives: 


Var{zrf}  =  /2</Var{i?2x} 


Equation  (3.21)  gives: 


Var{i?2X2}  =  Q2  exp(cr) 


m  + 1 


/« 


expcr2-! 


=  exp(<j2) 


m  + 1 


m 


expcr2-! 


since  Q2  =  1 . 
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Equating  Equations  (3.17)  and  (3.20)  and  taking  into  account  Equation  (3.21) 
gives  the  following: 


f2d  exp  a2 
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f  ,2  d) 
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l  m) 
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Applying  (3.51)  in  (3.52)  gives: 
m  + 1 
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Probability  of  bit  error 


APPENDIX  III-D.  MODELED  AND  SIMULATED  PROBABILITY 
OF  BIT  ERROR  FOR  THE  NAKAGAMI-LOGNORMAL  CHANNEL 
USING  60°  ANTENNA  SECTORING 


Figure  3.58.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  0.5) 
and  Lognonnal  Shadowing  ( adB  =  4)  Using  Six  Sectors  and  a  Rate  1/2  Convolutional 

Encoder  with  v  =  8. 
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Probability  of  bit  error 


Figure  3.59.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.75) 
and  Lognonnal  Shadowing  ( crdB  =  4)  Using  Six  Sectors  and  a  Rate  1/2  Convolutional 

Encoder  with  v  =  8. 


84 


Probability  of  bit  error 


Figure  3.60.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  1)  and 
Lognormal  Shadowing  ( crdB  =  4)  Using  Six  Sectors  and  a  Rate  1/2  Convolutional 

Encoder  with  v  =  8. 
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Probability  of  bit  error 


Figure  3.61.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  1.5) 
and  Lognonnal  Shadowing  ( adB  =  4)  Using  Six  Sectors  and  a  Rate  1/2  Convolutional 

Encoder  with  v  =  8. 
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Probability  of  bit  error 


Figure  3.62.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  2)  and 
Lognormal  Shadowing  ( crdB  =  4)  Using  Six  Sectors  and  a  Rate  1/2  Convolutional 

Encoder  with  v  =  8. 
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Probability  of  bit  error 


APPENDIX  III-E.  PROBABILITY  OF  BIT  ERROR  FOR  THE 
NAKAGAMI-LOGNORMAL  CHANNEL  KEEPING  THE  NUMBER 
OF  USERS  FIXED  AND  VARYING  THE  NAKAGAMI  PARAMETER 


SNR  per  bit  (dB) 

Figure  3.63.  Probability  of  Bit  Error  for  DS-CDMA  with  FEC  in  Various  Fading 
Conditions  with  40  and  90  Users  Per  Cell,  60°  Sectoring,  Lognormal  Shadowing 
( (JdB  =  4)  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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Probability  of  bit  error 


Figure  3.64.  Probability  of  Bit  Error  for  DS-CDMA  with  FEC  in  Various  Fading 
Conditions  with  120  Users  Per  Cell,  60°  Sectoring,  Lognormal  Shadowing  (adH  =  4) 
Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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Figure  3.65.  Probability  of  Bit  Error  for  DS-CDMA  Using  Sectoring  for  {crdB  =  4) ,  m 
=  0.5  and  10  Users  Per  Adjacent  Cell  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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IV.  EFFECT  OF  USER  DISTRIBUTION  AND  POWER  CONTROL 
ON  THE  FORWARD  CHANNEL  MODEL 


In  Chapter  III,  the  performance  of  the  DS-CDMA  cellular  system  in  a  Nakagami 
fading  environment  including  Lognormal  Shadowing  for  different  values  of  the 
Nakagami-/;/  parameter  was  analyzed.  This  analysis  assumes  that  the  intended  user  was 
in  the  least  optimum  position  within  the  cell,  which  can  be  any  comer  of  the  hexagon  as 
shown  in  Figure  3.1.  The  system  perfonnance  can  be  improved  by  assuming  that  the 
intended  user  can  be  anywhere  in  the  cell,  and  whose  position  can  be  represented  by  a 
random  variable  from  a  specific  probability  density  function.  A  revised  version  of  bit 
error  probability  can  be  defined  by  incorporating  this  probability  density  function  in  the 
initial  formula  for  the  bit  error  probability.  This  procedure  will  be  examined  in  Section 
IV-A.  Additionally,  the  analysis  in  Chapter  III  assumed  that  the  base  station  was 
transmitting  a  fixed  equal  power  to  all  users.  In  Section  IV-B,  the  base  station  transmits 
more  power  to  users  who  need  it  and  less  power  to  users  who  do  not  need  it.  System 
perfonnance  is  thus  enhanced  because  the  co-channel  interference  is  reduced.  The 
analysis  is  similar  to  that  given  in  [1],  For  simplicity  reasons,  the  hexagonal  cells  will  be 
replaced  by  overlapping  circular  cells  as  depicted  in  Figure  4.1.  Overlapping  is  used  for 
soft-handoff  reasons  of  mobile  users  between  cells.  Referring  to  Figure  4.1,  the  center 
base  station  is  located  at  the  origin  of  the  polar  coordinate  system,  with  the  intended 
mobile  user’s  position  in  the  center  cell  to  be  defined  by  ( r,0 )  where  0  <  r  <1  and 
-it  <9  <n . 
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A 


l 

* 


Figure  4.1.  Circular  Seven-Cell  Cluster  from  [1], 

Based  on  the  previous  analysis  in  Chapter  III,  it  is  known  that  the  distance 
between  the  base  station  in  the  center  cell  to  each  of  the  adjacent  base  stations  is  -J3 
since  d  =  1  is  normalized.  The  distance  Dj  from  each  adjacent  base  station  to  the  mobile 
user  in  the  center  cell  as  a  function  of  r  and  9  is  given  in  [1]  as  follows: 


-Jr 2  -  2-j3r  cos  9  +  3, 

i  =  1 

-Jr 2  -  V3 r  cos  9  -  3 r  sin  0  +  3, 

i  =  2 

yjr2  +  \[3r  cos  9  -  3r  sin  9  +  3, 

i=3 

-Jr2  +  2 yf3r  cos  0  +  3, 

i  =  4 

-Jr2  +  -J3r  cos  9  +  3r  sin  9  +  3, 

i  =  5 

-Jr2  - -J3r  cos  0  +  3r  sin  9  +  3, 

i  -  6 

These  distances  are  used  in  order  to  calculate  the  path  loss  for  the  intended 
information  signal  and  the  co-channel  interference.  Since  these  distances  D.  are 
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randomized,  it  is  too  difficult  and  complicated  to  use  the  Hata  model  in  order  to  predict 
the  path  loss.  A  more  convenient  and  practical  fonnula  to  compute  the  path  loss  is  as 
follows: 


Ln{d)' 


V^o  j 


(4.1) 


where  n  =  4  is  the  path  loss  exponent. 

Following  the  notation  as  in  [1],  and  assuming  that  the  intended  mobile  user’s 
position  is  fixed  and  exactly  defined  using  the  parameters  r  and  0  as  rf  and  0 1 

respectively,  (3.3)  can  be  revised  to 
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it  M  f  Di  2  Ec 


(4.2) 


where  Ec  is  given  as 


E  =  Wcc  _  fPtTcc 


L,\rf) 


(4.3) 


The  use  of  a  instead  of  a  in  (3.4)  gives  similar  results  for  bit  error  probability  as 
depicted  in  [1]  in  the  case  of  equal  transmitting  power  for  all  user  channels. 


USER  DISTRIBUTION  IN  THE  CELL 


As  mentioned  at  the  beginning  of  Chapter  IV,  it  is  possible  to  enhance  the 
performance  of  the  DS-CDMA  cellular  system,  assuming  that  the  users  are  physically 
distributed  in  the  cell  following  a  specific  probability  distribution  function  rather  than  the 
intended  mobile  user  being  placed  in  the  worst  case  position  in  the  cell.  In  the  case  that 
the  user  density  is  uniform  and  assuming  that  the  probability  density  function  pe{9 )  for 
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the  angular  positions  of  a  user  is  uniformly  distributed,  the  first  event  error  probability  is 
developed  in  [1]  and  is  given  as  follows: 


~  rn  pi  v 

W)=f  f  f  -Q 

J  —n  Jo  J  -oo 


exp 


f  2  A 
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6  K<- 1 


3  N 


I  4  r 


-  +  ■ 
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i=l  j= o 


fx  Di  2  E/l 


PzLi  ( zd)dzddrdO  (4.4) 


where  zd  is  the  value  of  the  random  variable  Zd  which  is  the  sum  of  d  Nakagami- 
square-Lognormal  random  variables. 

Substituting  (4.4)  in  (3.5),  an  upper  bound  on  the  bit  error  probability  for  the 
N a k a g a m i - Lo g n o rm a  1  fading  channel  is  obtained.  The  integral  in  (4.4)  was  simulated 
using  95,000  Monte  Carlo  trials  for  different  values  of  the  Nakagami-m  parameter  with 
60°  antenna  sectoring.  Figures  4.2  through  4.6  indicate  how  the  system  performs 
compared  to  the  worst  case  user  position  analyzed  in  Chapter  III  in  different  fading 
conditions  with  a  Lognonnal  Shadowing  value  of  crdB  =  7  . 
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Probability  of  bit  error 


Figure  4.2.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.5) 
and  Lognonnal  Shadowing  ( adB  =  7  )  Applying  Linear  User  Distribution  with  60° 
Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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Probability  of  bit  error 


Worst-Case  SNR  per  bit  (dB) 

Figure  4.3.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  0.75) 
and  Lognormal  Shadowing  ( aiB  =  7  )  Applying  Linear  User  Distribution  with  60° 
Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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Probability  of  bit  error 


Figure  4.4.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =1)  and 
Lognormal  Shadowing  ( (JdB  =  7  )  Applying  Linear  User  Distribution  with  60°  Sectoring 
and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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Probability  of  bit  error 


Figure  4.5.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  1.5) 
and  Lognonnal  Shadowing  ( adB  =  7  )  Applying  Linear  User  Distribution  with  60° 
Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 
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Figure  4.6.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  2)  and 
Lognormal  Shadowing  ( adB  =  7  )  Applying  Linear  User  Distribution  with  60°  Sectoring 
and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with  v  =  8. 


Observing  Figures  4.2  through  4.6,  it  is  obvious  that  the  system  performs  better  (3 
to  4  dB)  in  various  fading  environments  when  the  linear  user  distribution  compared  to  the 
worst-case  position  of  the  user  is  incorporated. 


B.  POWER  CONTROL  ON  THE  FORWARD  CHANNEL 


In  Section  A  the  user  distribution  in  the  performance  analysis  was  incorporated 
and  better  results  were  obtained  compared  to  those  in  the  previous  chapter.  The  system 
performance  can  be  enhanced  even  more  if,  in  addition  to  user  distribution,  the  power 
received  by  each  user  is  adjusted  to  be  equal.  Working  in  this  manner,  the  transmitting 
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power  to  non-intended  users  is  reduced,  which  leads  to  a  reduction  of  co-channel 
interference  from  users  located  in  other  cells. 


In  previous  chapters,  the  baseline  transmit  power  was  defined  as  Pt .  Assume 

that: 


P  = 


a 


(4.5) 


is  the  target  level  power  which  each  user  should  receive  and  a  is  an  attenuation  factor. 
The  values  in  (4.5)  are  constant  for  all  users  in  all  cells.  Each  mobile  user,  when 
receiving  the  information  signal,  measures  the  actual  power  received  and  reports  it  back 
to  the  base  station.  From  that  point  on,  the  received  power  for  each  mobile  user  is  not  a 
random  variable  as  mentioned  in  (2.12).  It  can  stated  that  it  is  a  realization  of  the  random 
variable  as  described  in  [1]: 


-  -  A  _  _  M  -  f„“P 

L„  ( dk )  xk  Ln  ( dk )  xk  dk  xk 


(4.6) 


where  dk  is  the  user’s  actual  distance  and  xk  is  the  shadowing  experienced  by  the 

specific  user.  In  order  for  each  user  to  receive  the  same  target  power  level  P,  the  base 
station  has  to  adjust  the  power  factor  as  follows: 

/»=^L  (4.7) 

a 

Substituting  (4.7)  into  (4.6)  implies: 

Pk=P  (4.8) 

Adjusting  the  power  factor  to  achieve  the  target  power  level  P  for  all  users  makes  it 
possible  to  overcome  the  Lognormal  Shadowing  effect  on  the  information  signal. 

Thus,  [1]  gives  the  first  event  error  probability  P2(d)  conditioned  or  rn  r,  6  as 
follows: 
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(4.9) 


h(d) 


rhr,g 


where  rf  is  the  value  of  the  Nakagami-square  random  variable  and  /;,  as  in  [1]  is  the 

random  variable  representing  the  distance  between  mobile  user  j  in  adjacent  cell  i,  with 
the  same  pdf  as  the  random  variable  r.  Recalling  the  probability  density  function  of  the 
sum  of  d  Nakagami-/;;  square  random  variables  from  (3.9)  and  the  fact  that 

pR  g  ( r,9 )  =  — ,  the  dependency  on  r ; ,  r,  6  can  be  removed  as  follows: 


Aw=f  r  f  e 


exP{i^ly  ^  Efe'}  ,  N0 
3 N  if  p0  D”  2 Ec 


Pu(u) 


(4.10) 


x  pRg(r,d)drdddu 

d 

where  u  =  ^  if  is  the  value  of  the  sum  of  d  Nakagami-square  random  variables.  A 

/= i 

closed  fonn  solution  for P2  (d)  is  given  as  follows: 


f1  r 

f  Q  1 

J0 

v  lam  +  0/ 

T(md  +  ^) 


2yfnT(jnd  +  1) 


^  \^md 

lam  +  Q 
2am 


(4.11) 


x 


i- 

k= 0 


(md  +  -)* 


lam 
lam  +  Q 


(md  + 1)/; 


-drdd 
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where  a- 


exp{i  2(72dB} 
3 N 


6  K- 1 

I  i 


i= 1  y=0 


Efo'l ,  k 

d;  2  Er 


Analytical  computation  of  the  closed  fonn  solution  can  be  found  in  Appendix  IV-A. 
Substituting  (4.10)  or  (4.11)  in  (3.5)  it  is  possible  to  compute  an  upper  bound  on  the 
probability  of  bit  error.  The  integral  in  (4.10)  was  simulated  using  95,000  Monte  Carlo 
trials  for  different  values  of  the  Nakagami-/?z  parameter  with  60°  antenna  sectoring.  The 
average  received  SNR  per  bit  for  the  Nakagami-Lognormal  channel  with  perfect  power 
channel  is  from  [1]  as  follows: 


n  =  e 


(4.12) 


As  shown  in  Figure  4.7,  in  the  case  of  ten  users  per  cell,  the  power-controlled 
system  works  better  than  the  system  with  fixed  power  (~6  dB).  However,  in  the  case  of 
thirty  users  per  cell,  the  power  controlled  system  works  better  than  the  fixed  power 
system  for  yb  <  1 7  dB  which  is  the  point  where  the  interference  floor  seems  to  be 

developed.  This  phenomenon  is  explained  analytically  in  [1]  and  in  a  few  words  means 
that  by  increasing  the  received  power  for  all  users  in  order  to  overcome  the  Lognormal 
Shadowing  effect,  the  co-channel  interference  is  simultaneously  increased.  Furthermore, 
in  the  case  of  m  =  0.5,  the  system  seems  to  be  able  to  accommodate  up  to  30  users  with  a 
bit  error  probability  of  ~  10“4  in  the  region  10  dB  <  fh  <  17  dB  . 
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Figure  4.7.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.5) 
and  Lognormal  Shadowing  (odB  =7)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 


Similar  results  are  obtained  when  changing  from  m  =  0.5  to  m  =  1.5  as  seen  in 
Figure  4.8.  In  that  case,  the  system  accommodates  up  to  60  users  per  cell  with  a  bit  error 
probability  of  ~  10“4  when  yh  <  16  dB . 
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Figure  4.8.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  1.5) 
and  Lognormal  Shadowing  (adfj  =7)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 


When  the  system  operates  in  more  lightly  shadowing  enviromnents  ((Jdfj  <  6) ,  the 

power-controlled  system  outperforms  the  fixed  power  system  as  depicted  in  Figure  4.9 
for  the  case  of  140  users  per  cell. 
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Figure  4.9.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  2)  and 
Lognormal  Shadowing  ( a dB  =4)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 


In  these  circumstances,  the  benefits  of  eliminating  the  Lognormal  Shadowing  of 
the  signal  overcomes  the  increase  to  co-channel  interference. 

Additionally,  system  designers  must  be  very  careful  when  the  system  operates  in 
heavily  Lognormal  Shadowing  environments  (odB  >  6)  in  order  for  the  power  controlled 
system  to  operate  out  of  the  interference-limited  region  even  when  using  60°  sectoring. 
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Figure  4. 10.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.5) 
and  Lognormal  Shadowing  (<7dfj  =7)  vs.  Users  per  Cell  Applying  Forward  Power 
Control  and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc=  1/2  and  v  =  8),  j^  =  15dB. 


In  the  case  of  m  =  0.5  and  a dB  -  7  >  6 ,  the  power  controlled  system  roughly 
accommodates  up  to  35  users  per  cell,  for  fh  =  15  dB  . 
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Figure  4.11.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  1)  and 
Lognormal  Shadowing  (<7dfj  =7)  vs.  Users  per  Cell  Applying  Forward  Power  Control 
and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc  =1/2  and  v  =  8),  y%15dB. 
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Figure  4.12.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  2)  and 
Lognormal  Shadowing  ((JdB  =  4)  vs.  Users  per  Cell  Applying  Forward  Power  Control 
and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc  =1/2  and  v  =  8),  y^  =  12dB. 


By  increasing  the  value  of  the  Nakagami-m  parameter  to  m=  1,  the  severity  of  the 
fading  is  decreased  and  thus  the  power-controlled  system  can  accommodate  up  to  64 
users  per  cell  for  fb  =  15  dB  as  shown  in  Figure  4.11.  When  m  =  2  and  crdB  =4  the 
power  controlled  system  outperforms  the  one  with  fixed  power. 

If  instead  of  six  sectors,  120°  sectoring  or  no  sectoring  is  used,  the  amount  of  co¬ 
channel  interference  will  be  increased  by  a  factor  of  two  and  six  respectively  which  will 
result  in  the  worst  performance  of  the  power  controlled  system.  In  that  case,  the  number 
of  users  per  cell  that  the  system  can  accommodate  in  the  no  interference-limited  region 
will  decrease. 

Figures  similar  to  the  aforementioned  figures  can  be  found  in  Appendix  IV-B  for 
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Lognormal  Shadowing  adB  -  4  and  adB  -  9 ,  which  reinforce  the  previous  comments. 
Especially  when  (JdB  =  9  and  in  =  0.5,  m  =  0.75,  m  =  1  the  fixed  power  system 

outperforms  the  one  with  power  control  as  depicted  in  Figures  4.21,  4.22  and  4.23  of 
Appendix  IV-B.  That  means  in  heavy  shadowing  and  severe  fading  environments  system 
designers  must  use  fixed  power  system  and  not  the  power  controlled  one. 

Accordingly,  it  can  be  stated  that  by  carefully  adding  power-control  to  the 
forward  channel,  it  is  possible  to  enhance  the  system  performance  operating  in  a  lightly 
shadowed  and  not  severe  fading  Nakagami-Lognormal  channel. 


C.  SUMMARY 


In  this  chapter,  the  system  perfonnance  was  enhanced  by  considering  that  the 
user’s  positions  in  the  cell  are  randomly  distributed  following  a  linear  distribution.  That 
assumption  is  more  realistic  compared  to  the  worst-case  position  of  the  intended  user 
located  at  the  corner  of  the  hexagon.  Given  the  revised  version  of  the  first  event  error 
probability  referring  to  the  more  realistic  case,  it  was  possible  to  compute  an  upper  bound 
on  the  probability  of  bit  error  by  simulating  the  integral  of  (4.4)  by  using  95,000  Monte 
Carlo  trials.  The  system  performance  was  improved  by  3  to  4  dB  for  udB  =  7  . 

Furthermore,  in  Section  IV-B,  a  power  control  technique  was  applied.  Thus,  all 
users  received  the  same  amount  of  power.  This  technique  overcomes  the  effects  of 
Fognormal  Shadowing  and  distance  dependent  path  loss.  Simulating  the  integral  of 
(4.10)  by  using  95,000  Monte  Carlo  trials,  an  upper  bound  on  the  probability  of  bit  error 
was  computed.  By  properly  using  the  power  control  technique,  it  was  possible  to  greatly 
increase  system  performance  and  user  capacity  in  the  forward  channel  operating  in  a 
Nakagami-Lognormal  Shadowing  environment.  As  mentioned  previously,  system 
designers  must  be  very  careful  when  incorporating  power  control  in  heavily  shadowed 
environments  in  order  for  the  system  to  operate  in  the  no  interference  limited  region. 
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APPENDIX  IV-A.  ANALYTICAL  COMPUTATION  OF  A  CLOSED 
FORM  SOLUTION  FOR  THE  FIRST  EVENT  ERROR 
PROBABILITY  IN  THE  CASE  OF  FORWARD  POWER  CONTROL 


Equation  (4.10)  gives  that: 


/u-o-r  r  f  g 


x  pRe(r,d)drdddu 


i 


exp{iV;*}£  g  E{r;}  |  jYo 


3  N 


i= 1  y=0 


D"  2E 


Pu{u) 


r  -  exp{i2cr^j  6  E{r"j  yy 

Substituting  pRe(r,0)  =  —  and  a- - ^  ^  - ^-  +  — —  we  get  ioxP2{d) 

Tt  3  N  i=1  y-()  E);  2^e 


the  following: 


w=i:  i;  i:  e 


n 


pu  (w)  dr dO du 


(4.13) 


Equation  (3.9)  gives: 


Pu{u) 


f  m  ^ 


md 


md-l 


T(md) 


exp 


f  nzw  ^ 

~o~ 


Applying  (4. 14)  to  (4.13)  we  get  for P2{d)  the  following: 

/V  C°°  1*1 

P2(d)=  f  J  f  Q 

J- oo  J-7T  Jo 


/  md-l 

m  u 

J 


r(md) 


r  mu  ^ 
x  exp - 


drdddu 


(4.14) 


(4.15) 


From  [9]  we  have  the  following  equality: 
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(4.16) 


n/4 


—y — -  f exp(-a?t)t/!  'Q ( \[ct) dt  =  - 1= - — - 

r(/?)J0  3;  ’  v+w  i4nr(p+\)(\  +  i//) 


x2Fi 


f  1  i  N 

hP+— ;/?+!;  - — 

2  l  +  ^J 


where  2Fl  is  the  Gaussian  hypergeometric  function  defined  as  follows: 


(  ,  6  -M** 

2  V  2  "  J  6  (/?+!)* 


(4.17) 


2  1  ^ 
iff  =  — ,  z  = - and 


1  +  ^ 


M*  = 


F(g  +  A:) 


r  U) 

g(g  +  l)...(g  +  A:-l) 


(&)o  = 1 

Comparing  now  the  variables  of  Equations  (4.15)  and  (4.16)  and  noting  the  random 

YYl  1 

variable  U  takes  only  positive  values  we  get  that:  a:,  = — ,/3  =  md,  t  =  w,  c  =  —  and 

£2  a 

Q 

\jf  =  — — .  Applying  these  last  five  equations  to  (4.14)  and  making  use  of  (4.15)  we  get 
lam 

for P2(d)  the  following: 


F  f1  r 

f  Q  1 

Jo  n 

v  lam  +  £iy 

r  (md  +  ^) 


x2Fx 


,  ,  1  ,  .  2  am 

1,  md  +  — ;  md  + 1;  — — - 

2  lam  +  £2 


l\fnT{md  + 1) 
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f  x^md 

lam  +  Q 
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(4.18) 


dr dO 
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Substituting  1F]  by  its  definition,  we  finally  get  forP2(J)the  following  close  form 
solution: 


(4.19) 


(md  + \)k 


*1 
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APPENDIX  IV-B.  PROBABILITY  OF  BIT  ERROR  FOR  THE 
NAKAGAMI-LOGNORMAL  CHANNEL  USING  FORWARD 
CHANNEL  POWER  CONTROL 


Figure  4.13.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  0.5) 
and  Lognormal  Shadowing  (crdB  =  4)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 
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Figure  4.14.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.75) 
and  Lognormal  Shadowing  (crdB  =  4)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 
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Figure  4.15.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  1)  and 
Lognormal  Shadowing  (c rdB  =4)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 
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Figure  4.16.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  1 .5) 
and  Lognormal  Shadowing  (crdB  =  4)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 
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Figure  4.17.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  2)  and 
Lognormal  Shadowing  (c rdB  =4)  Applying  Forward  Power  Control  and  Linear  User 

Distribution  with  60°  Sectoring  and  FEC  Using  a  Rate  1/2  Convolutional  Encoder  with 

v  =  8. 
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Figure  4.18.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  3)  and 
Lognormal  Shadowing  [adB  =  4)  vs.  Users  per  Cell  Applying  Forward  Power  Control 
and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc  =1/2  and  v  =  8),  ^  =  12dB. 
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Figure  4.19.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  4)  and 
Lognormal  Shadowing  [adB  =  4)  vs.  Users  per  Cell  Applying  Forward  Power  Control 
and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc  =1/2  and  v  =  8),  ^  =  12dB. 
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Figure  4.20.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.75) 
and  Lognormal  Shadowing  {odft  =  7)  vs.  Users  per  Cell  Applying  Forward  Power 
Control  and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc  =1/2  and  v  =  8),  ^  =  15  dB. 
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Figure  4.21.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  0.5) 
and  Lognormal  Shadowing  [adH  =  9)  vs.  Users  per  Cell  Applying  Forward  Power 
Control  and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(f?cc  =1/2  and  v  =  8),  y~b  =  20  dB  . 
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Figure  4.22.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  ( m  =  0.75) 
and  Lognormal  Shadowing  (rx^  =  9)  vs.  Users  per  Cell  Applying  Forward  Power 
Control  and  Linear  User  Distribution  using  with  60°  Sectoring  and  FEC 
[Rcc  =1/2  and  v  =  8),  y~b  =  20  dB  . 
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Figure  4.23.  Probability  of  Bit  Error  for  DS-CDMA  with  Nakagami  Fading  (m  =  1)  and 
Lognormal  Shadowing  (adB  =  9)  vs.  Users  per  Cell  Applying  Forward  Power  Control 
and  Linear  User  Distribution  with  60°  Sectoring  and  FEC 
(Rcc=  1/2  and  v  =  8),  y~b  =  20  dB  . 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


In  this  thesis,  it  was  possible  to  look  at  the  performance  on  the  forward  channel  of 
a  six-sectored  DS-CDMA  cellular  system  with  FEC  operating  in  a  Nakagami  fading  and 
lognormal  shadowing  environment.  Attempts  were  made  to  enhance  the  perfonnance  of 
the  system  by  using  different  techniques  such  as  linear  user  distribution  and  power 
control. 

A.  CONCLUSIONS 


In  Chapter  II,  a  forward  channel  for  the  DS-CDMA  cellular  system  was  created. 
By  taking  into  account  all  the  appropriate  losses  predicted  by  the  extended  Hata  model 
and  using  Nakagami  fading  in  order  to  include  small-scale  propagation  effects,  the  model 
for  the  transmitted  signal  for  the  forward  channel  was  built.  Additionally,  Lognonnal 
Shadowing  was  incorporated.  The  result  for  a  model  that  included  all  the  aforementioned 
characteristics  of  the  propagation  channel  was  a  closed  fonnula,  which  describes  the 
received  signal  to  a  mobile  user’s  receiver. 

In  Chapter  III,  given  the  expression  for  the  upper  bound  on  the  probability  of  bit 
error  by  using  Forward  Error  Correction  (FEC)  in  the  fonn  of  convolutional  encoding 
with  soft-decision  decoding,  the  performance  of  the  cellular  system  was  analyzed  in  two 
ways.  A  statistical  model  Zd  was  created  in  order  to  approximate  the  sum  of  d 

multiplicative  Nakagami-square -Lognonnal  random  variables  as  a  multiplicative 
Nakagami-square-Lognormal  random  variable.  By  using  this  statistical  model,  it  was 
possible  to  compute  an  upper  bound  of  the  probability  of  bit  error  in  different  fading 
environments  varying  the  Nakagami -m  parameter  from  0.5  to  2.  These  results  were 
compared  with  the  Monte  Carlo  simulated  results.  This  model  was  accurate  in  the  SNR 
per  bit  range  from  yh  =10  to  15  dB  for  all  values  of  in.  Specifically,  the  model  was 

tested  for  m  =  0.5,  m  =  0.75,  m  =  1,  m  =  1.5  and  m  =  2.  These  specific  values  of  m  vary 
the  severity  of  the  fading  channel  from  one-sided  Gaussian  fading  (in  =  0.5)  to  Ricean 

fading  (in  =  2).  The  analysis  in  this  chapter  took  into  account  the  intended  user  was  at  the 
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worst-case  position  in  the  corner  of  the  center  cell.  The  model  is  a  very  helpful  tool  for 
system  designers,  which  can  be  used  quickly  and  effectively  for  analyzing  the 
performance  of  DS-CDMA  cellular  systems  operating  in  a  Nakagami-Lognonnal 
channel. 

Also  in  this  chapter,  the  decrease  in  the  system’s  perfonnance  was  mentioned 
when  using  three  sectors  instead  of  six,  or  the  no  sectoring  case.  Furthermore,  the 
system’s  performance  was  examined  with  different  values  of  users  per  adjacent  cell  and 
for  Lognormal  Shadowing  7  dB  and  4  dB.  The  enhancement  in  the  system’s 
performance  from  7  dB  to  4  dB  was  very  obvious.  (Lightly  shadowing  environment.) 

In  Chapter  IV,  the  performance  analysis  in  Chapter  III  was  modified  by 
incorporating  user  distribution  within  the  cell  and  assuming  that  the  user’s  position  is 
random  following  a  specified  distribution.  Using  a  modified  formula  for  the  upper  bound 
on  the  probability  of  bit  error  for  the  DS-CDMA  cellular  system  with  FEC  in  Nakagami- 
Lognonnal  channel  in  order  to  incorporate  the  user’s  distribution,  the  performance  of  the 
system  was  analyzed  with  the  new  assumption  included.  It  is  worth  mentioning  that  in 
this  specific  case,  the  system  performance  was  3  to  4  dB  better  than  predicted  in  Chapter 
III,  the  worst-case  scenario,  using  the  results  predicted  by  the  Monte  Carlo  simulations  in 
both  cases.  Additionally,  power  control  was  implemented  in  the  forward  traffic  channel 
in  a  way  that  each  user  received  the  same  amount  of  power.  That  fact  overcame  the 
effects  of  Lognonnal  Shadowing  and  distance  dependent  path  loss.  Once  again  using 
Monte  Carlo  trials,  it  was  discovered  that  by  properly  using  the  power  control  technique, 
the  system  performance  was  greatly  improved.  On  the  other  hand  however,  system 
designers  must  be  very  careful  when  incorporating  power  control  in  the  forward  channel 
in  order  for  the  system  to  operate  in  the  no  interference  limited  region  especially  in  heavy 
shadowing  environments. 


B.  FUTURE  WORK 


It  is  possible  to  further  analyze  the  perfonnance  of  DS-CDMA  cellular  systems 
operating  in  a  Nakagami  fading  and  lognormal  shadowing  channel.  For  example,  FEC 
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can  be  used  with  Serially  Concatenated  Convolutional  Codes  (SCCC),  which  is  a  more 
sophisticated  type  of  coding  in  order  to  improve  system  performance. 

Furthermore,  an  error  at  the  pilot  tone  recovery  can  be  introduced  which  may 
yield  more  realistic  results.  Additionally,  in  Chapter  IV,  different  types  of  user 
distributions  can  be  used  which  may  produce  better  results. 

Finally,  a  similar  statistical  model  can  be  applied  to  the  reverse  channel  of  a  DS- 
CDMA  cellular  system  in  a  Nakagami  fading  and  Lognormal  Shadowing  environment. 
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